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Abstract 

We reconsider the effective Lagrangian that describes a light Higgs-like boson and bet- 
ter clarify a few issues which were not exhaustively addressed in the previous literature. 
In particular we highlight the strategy to determine whether the dynamics responsi- 
ble for the electroweak symmetry breaking is weakly or strongly interacting. We also 
discuss how the effective Lagrangian can be implemented into automatic tools for the 
calculation of Higgs decay rates and production cross sections. Finally, we present 
eHDECAY, a modified version of the program HDECAY which includes the full list of 
leading bosonic operators of the effective Lagrangian and implements two benchmark 
composite Higgs models. 



1 Introduction 



The exploration of the weak scale has marked an important step forward with the discovery 
by the ATLAS [l] and CMS [2] collaborations of a boson with mass rrih — 125 GeV, whose 
production cross section and decay rates are compatible with those predicted for the Higgs 
boson of the Standard Model (SM). At the same time, no hint of the existence of additional 
new particles has emerged yet, which might shed light on the origin of the electroweak sym- 
metry breaking (EWSB). One is thus faced with the problem of which is the best strategy to 
describe the properties and investigate the nature of the new boson h, beyond the framework 
of the Standard Model. In absence of a direct observation of new states, our ignorance of the 
EWSB sector can be parametrized in terms of an effective Lagrangian for the light boson. 
Such an effective description is valid as long as New Physics (NP) states appear at a scale 
M 3> rrih, and is based on an expansion in the number of fields and derivatives [3||. The 
detailed form of the effective Lagrangian depends on which assumptions are made. Consider- 
ing that the observation made by the LHC experiments is in remarkable agreement with the 
SM prediction, although within the current limited experimental precision, it is reasonable 
to assume that h is a CP-even scalar that forms an SU(2)l doublet together with the lon- 
gitudinal polarizations of the W and Z, so that the SU(2)l X U(l)y electroweak symmetry 
is linearly realized at high energies. Under these assumptions the effective Lagrangian can 
be expanded into a sum of operators with increasing dimensionality, where the leading NP 
effects are given by dimension-6 operators. 

We reconsider such an effective Lagrangian in Section |2j By means of a naive power 
counting we estimate the coefficients of the various operators and review the most impor- 
tant bounds set on them by present experimental results on electroweak (EW) and flavor 
observables. Focusing on Higgs physics, we then discuss in Section [3] the relative effect of 
the various operators on physical observables. Such an analysis, first proposed in Ref. [2], 
will allow us to identify which operators can probe the Higgs coupling strength to the new 
states and which instead are sensitive only to the mass scale M. This is of key impor- 
tance to distinguish between weakly-coupled UV completions of the Standard Model, like 
Supersymmetric (SUSY) theories, and theories where the EW symmetry is broken by a new 
strongly-interacting dynamics which forms the Higgs boson as a bound state [5]-[7j|4] . These 
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are the two most compelling scenarios put forward to solve the hierarchy problem of the 
Standard Model. We conclude the section by discussing how the assumption of a Higgs 
doublet and linearly-realized SU(2)l x U(l)y can be relaxed. We illustrate the non-linear 
effective Lagrangian valid for the case of a generic CP-even scalar h and discuss the impli- 
cations of custodial invariance. Section [4] is devoted to clarify a few issues related to the use 
of the effective Lagrangian beyond the tree level. We present our concluding discussion in 
Section |5j In the Appendices [A|fC] we collect useful formulas and give further details on the 
construction of the effective Lagrangian. The details of how we derived the bounds on the 
dimension-6 operators are reported in Appendix D As an illustration of our analysis and to 
better demonstrate how the effective Lagrangian can be implemented into automatic tools 
for the computation of physical quantities like Higgs production cross sections and decay 
rates, in Appendix [E] we present eHDECAYQ a modified version of the program HDECAY which 
includes the full list of leading bosonic operators. 

2 Effective Lagrangian for a light Higgs doublet 

The most general SU (3)c x SU (2)lXU (l)y-invariant Lagrangian for a weak doublet H at the 
level of dimension-6 operators was first classified in a systematic way in Refs. [8] . Subsequent 
analyses [9,10 pointed out the presence of some redundant operators, and a minimal and 



complete list of operators was finally provided in Ref . [Tl] . As recently discussed in Ref . (2] , 
a convenient basis of operators relevant for Higgs physics is the following: 

C = C SM + = C sm + ACsiLH + A£ Fl + ACp 2 (2.1) 



1 eHDECAY is available at the following URL: http://www-itp.particle.uni-karlsruhe.de/~maggie/ 
eHDECAY/ 
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with 

AC SILH = * d"(H'H) 3^H) + * (fltfi»fl) (jfttf ^ _ 5^ (fltfi)' 

+ f % V» H*Hq L H'u R +%y d H*Hq L Hd R + % y, H*H L L Hl R + h.c) 

+ 'ti, {^' 9h ) v™** + ^ 9h ) ^ (2.2) 

+ {D"H)^a\D v H)Wi u + (D^H)\D U H)B, U 

- /2 - 2 

+ ^f- rfHB^B^ + ^ H^HG a ^G a ^ u , 

272,^/ 



+ (w^*) (# ct ^#) + h.c)) (2.3) 

+ ^ {L L >fL L ) (H^,H) + ^ (L i7 VL L ) (HW^.H) 

A n CuB 9 - jj C nv r> i ^uW 9 - i TTC av rtri i ^uG9S - TTC au\a ria 

AC F2 = — 2~ l/« <?l# M fl H 2~ f« ?lct H a* u R H ^— y u q L H cx M A u R 

m w m w m w 

+ CjE f Vd q L Ha^d R + ^ y d q L a*Ha^d R + y d q L Ha^\ a d R G a 

+ C -^fyi L L Han R B^ + ^ Vl L L a l Ha^l R W% + h.c. 
m w m w 

(2.4) 

The SM Lagrangian Csm and our convention for the covariant derivatives and the gauge 
field strengths are reported for completeness in Appendix |A} In particular, A is the Higgs 
quartic coupling and the weak scale is defined to be 

' 246 GeV. (2.5) 
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By iH^9H we denote the Hermitian derivative iH\D^H)-i{D^H)^ H , and = i^, 7"]/2. 



The Yukawa couplings y u ,d,i and the Wilson coefficients q in Eq. (2.3) are matrices in flavor 



space, and a sum over flavors has been left understood. As specified in Eq. (2.1), we will 
denote as Oi the dimension-6 operator whose coefficient is proportional to Cj. 

Our higher-dimensional Lagrangian, which is supposed to capture the leading New Physics 
effects, counts 12 (ACsilh) + 8 (A£pJ + 8 (ACp 2 ) — 28 operators. Five extra bosonic oper- 
ators, 



C3W 9 Ajk-tjii vrxri pjirk a ^3G9s fabc/^tav /~ib p/^ca 

— —t J WfWp' 1 , —2—/ W ^/G/ 

IIL W L w 



m 



{D p W p ») 1 , 



C2B 

ml 



(d»B, u ) {d p B 



pv\ 



P ' 



(2.6) 



which affect the gauge-boson propagators and self-interactions but with no effect on Higgs 
physics, should also be added to complete the operator basis, as well as 22 four-Fermi baryon- 
number-conserving operators. ^ A comparison with Ref. 



{D p G pu ) a 



11 



shows that two of our operators 



are actually redundant. As we shall explain in more detail in Section pi (see Eqs. (3.28) 



(3.29)), it is well known 10 13 that two particular linear combinations of the fermionic 
operators in AC Fl are equivalent to pure oblique corrections parametrized by the operators 
Or-, Ow and Ob'- 



0^ = J2 Y ^°h^-O t ,O b and 0' Hq + 0' HL 



O 



W 



(2.7) 



where the sum runs over all fermion representations, ip = qL,UR,dR, L L ,l R , whose hyper- 
charge has been denoted as Y^. These two linear combinations have then to be excluded 
from ACf 1} and we end up with exactly 53 linearly- independent operators as in Ref. [IT]. ^ 
Any other dimension-six operator can be obtained from these 53 operators by using the 



2 Notice that the last three operators in Eq. (2.6 1 can be rewritten in favor of three additional independent 
four- Fermi operators, as in the basis of Ref. 



11 



. The coefficients C2W1 C 2B contribute respectively to the W 



12 



and Y parameters defined in Ref. 

3 For completeness we collect in Appendix [c] also the extra 6 bosonic operators of dimension-six that are 

CP-odd. 
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equations of motion, or equivalently by performing appropriate field redefinitions. 

Even though our basis (2.2)-(2.4) is equivalent to the one proposed in Ref. |11|, we 
advocate that it is more appropriate for Higgs physics for at least three reasons [4 : i) Generic 
models of New Physics generate a contribution to the oblique S parameter 



14 



12 



at tree- 



level, which in the basis of Ref. [IT] would have to be encoded in the two fermionic operators 
H ^ and 0' H + 0' HL even in the absence of direct couplings between the SM fermions and 
the New Physics sector. There is an advantage in describing the oblique corrections in terms 



of the operators in (2.2) rather than in terms of the operators with fermionic currents, which 



generate vertex corrections and modify the Fermi constant, ii) The basis (2.1) isolates the 
contributions to the decays h — > 77 (from 7 ) and h —> (from 7 and Ohw ~ Ohb) that 
occur only at the radiative level in minimally coupled theories, in) Our basis of operators 
is more appropriate to establish the nature of the Higgs boson and determine the strength 
of its interactions. For example, as we shall explain momentarily, if the Higgs boson is a 
pseudo Nambu-Goldstone boson (pNGB) the coefficient of the operator 7 , hence the rate 



h — > 77, is suppressed, while in the basis of Ref. 11 this reflects into a cancellation in the 



linear combination 4c 7 + (cww — cwb) (cf. footnote |4j) . 

While a complete classification of the operators is essential, having a power counting to 
estimate their impact on physical observables, hence their relative importance, is equally 
crucial. In this sense a simple yet consequential observation was made in Ref. [2]: when 
expanding the effective Lagrangian in the number of fields and derivatives, any additional 
power of H is suppressed by a factor g*/M = 1/f, where < An denotes the coupling 
strength of the Higgs boson to New Physics states and M is their overall mass scale; any 
additional derivative instead costs a factor \jM. If the light Higgs boson is a composite 
state of the dynamics at the scale M, it is natural to expect c/* 3> 1, hence / <C M, which 
implies that operators with extra powers of H give the leading corrections to low-energy 



4 In particular, the following identities hold 

a 2 I 



2- H^H W; V W { » v = O ww =O w ~0 B + HB - O hw + - 1 

(2.8) 



ffVff wLb>» = o WB = o B -o HB --o 1 . 
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observables. On the other hand, in weakly-coupled completions of the Standard Model 
where g* ~ g, all operators with the same dimension can be equally important. A proper 
analysis of the experimental results through the language of the effective Lagrangian can thus 
give indication on whether the dynamics at the origin of electroweak symmetry breaking is 
weakly or strongly interacting. According to the power counting of Ref. [I], one naively 
estimates (ip = u, d, I, q, L) ^ 



v 2 \ ~fm 2 w \ _ _ _ _ „f m w 



- -/ n( X l v2 \ - n f^uX d v 2 \ _ _ _ ( m 2 w 

c ^,c^ ~ CI — — I , c Hud ~U\—^—\, C^ W ,C^B,C^ G ~ U\—^j 



(2.9) 



167T 2 / 2 , 

where denotes the coupling of a generic SM fermion ip to the new dynamics. It should be 
stressed that these estimates are valid at the UV scale M, at which the effective Lagrangian 
is matched onto explicit models. Renormalization effects between M and the E W scale mix 
operators with the same quantum numbers, and give in general subdominant corrections to 
the coefficients. We shall comment on these renormalization effects in Section HI Notice 
that the estimates of cw,b, cn^p, c' H ^ and ct apply when these coefficients are generated at 
tree-level. However, specific symmetry protections which might be at work in the UV theory, 
like for example i?-parity in SUSY theories, can force the leading corrections to arise at the 
1-loop level. 



Equation (2.9 ) suggests that in the case of a strongly-interacting light Higgs boson (SILH) 
the leading New Physics effects in Higgs observables are parametrized by the operators 
Oh,t,6,*i>, an d' if SM fermions couple strongly to the new dynamics, by the fermionic 



operators of Eq. (2.3) [41. Notice that, compared to the naive counting, CHW,HB, g ,-y are 
suppressed by an additional factor (g 2 /16ir 2 ). This is because the corresponding operators 
contribute to the coupling of on-shell photons and gluons to neutral particles and modify the 
gyromagnetic ratio of the W, and are thus generated only at the loop level in a minimally 



coupled theory. Similarly, the dipole operators of Eq. (2.4) are generated at the loop-level 



^Notice that our normalization differs from the one of Ref. ! |4 , and it is more convenient than the latter 



for a model-independent implementation of Eq. (2.2 1 in a computer program. The factor multiplying each 
operator in the effective Lagrangian has been conveniently defined such that the dependence on M and 
is fully encoded in the dimensionless coefficients 5$. 
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only, hence their estimates have an extra loop factor. 

A special and phenomenologically motivated case is represented by theories where the 
Higgs doublet is a composite Nambu-Goldstone (NG) boson of a spontaneously-broken sym- 
metry Q — > H of the strong dynamics |5]-[7j|4). For these models the scale / must be identified 
with the decay constant associated with the spontaneous breaking, and the naive estimate 
of the Wilson coefficients Cj is modified by the request of invariance under Q in the limit of 
vanishing explicit breaking. At the level of dimension-6 operators, 7 , O g , Oq, Ut d,i and the 



dipole operators of Eq. (2.4) violate the shift symmetry H l — >■ H' 1 + ( l (£ ! = const.) that 
is included as part of the Q/TL transformations. This means that they cannot be generated 
in absence of an explicit breaking of the global symmetry. It follows, in particular, that 
the naive estimates of the operators 7 and O g carry in this case an additional suppression 
factor [2], 

m 2 w \ 9% 



where g@ denotes any weak coupling that breaks the Goldstone symmetry (one of the SM 
weak couplings in minimal models, i.e. the SM gauge couplings or the Yukawa couplings). 
The operators 6 , O^, O^g, 0<pw, O^b have been defined so that their prefactor already 
includes one spurion coupling, precisely the Higgs quartic coupling A in 6 , and the Yukawa 
coupling in the other operators - indeed, both these couplings vanish for an exact NG 
boson. The estimates of the corresponding coefficients cq, c^, c^g, fyw, c^b are thus not 
modified. 



In writing Eq. (2.2 ) we have assumed that each of the operators U)( y is flavor-aligned with 
the corresponding fermion mass term, as required in order to avoid large Flavor-Changing 
Neutral Currents (FCNC) mediated by the tree-level exchange of the Higgs boson (see for 
example Ref. [15] for a natural way to obtain this alignment). This implies one coefficient 
for the up-type quarks (c u ), one for down-type quarks (q), and one for the charged leptons 
(q), i.e. the c U) d,i are proportional to the identity matrix in flavor space. 

2.1 Current bounds on flavor-preserving operators 

It is useful to review some of the most important constraints on the coefficients q that follow 
from current experimental results, such as electroweak precision tests, flavor data and low- 
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energy precision measurements. For simplicity, we focus on the bounds on flavor- conserving 
operators, keeping in mind that they can come also from flavor-changing processes. For a 
discussion of the bounds on flavor-violating operators see for example the recent review of 
Ref. [Ui as well as Ref. [l7|. 

Among the strongest bounds are those on operators that modify the vector-boson self- 
energies. The operator Ot, for example, violates the custodial symmetry [18] and contributes 



to the EW parameter e\ 19 . From the EW fit performed in Ref. 20 , it follows, with 95% 
probability, 

Aei = Ap = c T {m z ) , -1.5 x 1(T 3 < c T {m z ) < 2.2 x 1(T 3 . (2.11) 

Such a stringent bound can be more naturally satisfied by assuming that the dynamics at 
the scale M possesses an (at least approximate) SU(2)v custodial invariance. In this case 
ct{M) = 0, and a non-vanishing value will be generated through the renormalization-group 
(RG) flow of this Wilson coefficient down to m z in the presence of an explicit breaking of the 
custodial symmetry, as due for example to the Yukawa or hypercharge couplings. We will 
discuss these renormalization effects in more detail in Section HI Notice that all the other 
dimension-6 operators in the effective Lagrangian are (formally) custodially symmetric and 
their coefficients will not be suppressed at the scale M. ^\ The electroweak precision tests 
also imply a strong bound on Ow + Ob [4] , since this linear combination contributes to the 



parameter e 3 (19] . With 95% probability, one has 20 



Ae 3 = c w (m z ) + c B (m z ), -1.4 x < c w (m z ) + c B (m z ) < 1.9 x . (2.12) 



From the tree-level estimate of cw,b reported in Eq. (2.9), and assuming an approximate 
custodial invariance to suppress ct as explained above, it follows that Eqs. (2.11) and (2.12) 
set a lower bound M > a few TeV. This bound is quite robust and can be avoided only in 
weakly-coupled UV completions where an extra symmetry protection suppresses the leading 



6 More precisely, for all the other operators the only violation of the custodial symmetry comes from 
the explicit breaking due to the gauging of hypercharge. As such, this breaking is external to the EWSB 
dynamics, since it comes from the weak gauging of its global symmetries. Formal invariance of the operators 
can be restored by uplifting the hypercharge gauge field to a whole triplet of SU(2)r. The top Yukawa 
coupling is another source of explicit custodial breaking. 



8 



contribution to cw,b by an additional loop factor. Notable examples are SUSY theories with 
i?-parity. 



The fermionic operators in Eq. (2.3) are strongly constrained by Z-pole measurements, 
as they modify the couplings of the Z to quarks and leptons: 

9i*l> 2 T 3L~Q sm 2 8 w ' g Rl p 2 Q sin 2 ^ ' 
where T 3i and Q are respectively the SU{2) L and electric charges of the fermion ip, and 



^ = {L,q} is the SU(2)l doublet to which tpi, belongs. We used the results of Ref. 20 



to 



perform a fit on the coefficients cuip, ch^, c' h ^,. The details of our analysis can be found in 



Appendix D (see also Ref. [21] ). In the case of light quarks (u,d,s) we find the following 
bounds 

-0.03 < c Hq i < 0.02 , -0.002 < c Hql < 0.003 , 
-0.005 < c Hq2 < 0.003 , -0.003 < c' Hq2 < 0.005 , (2.14) 
-0.008 < c Hu < 0.02 , -0.03 < c Hd < 0.02 , -0.03 < c Hs < 0.02 , 
while a fit on leptons and heavy quarks (c, b) gives 

-0.004 < c HL + c HL < 0.002 , -0.003 < c HL - c HL < 0.0002 , -0.0007 < c m < 0.003 , 
-0.02 < c Hq2 + c Hq2 < 0.005 , -0.02 < c Hc < 0.03 , 
-0.003 < c Hq3 - c Hqs < 0.009 , -0.07 < c Hb < -0.005 . 

(2.15) 

All the above bounds have 95% probability and by the various coefficients we mean their 
values at the scale mz- The weakest constraint is that on the operator Onb, which modifies 
the coupling of 5r to the Z boson. The operator involving two right-handed top quarks, 
Ont, is unconstrained by EW data, but it is also not relevant for the Higgs decays and 
will be neglected in the following. The coefficient cntb is severely constrained by the 6—7-57 
rate. Indeed, the expansion of Outb around the vacuum contains a vertex of the type Wt^bR, 
which at 1-loop gives a chirally-enhanced contribution to the rate (see for example Ref. [22]). 
We find, with 95% probability: 

- 0.4 x 10" 3 < c Htb (m w ) < 1.3 x 10~ 3 . (2.16) 
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For a given (v/f), the above bounds set a limit on the couplings of the SM fermions to the 



new dynamics, see Eq. (2.9). Unless the scale of New Physics is very large, or some specific 



symmetry protection is at work in the UV theory (see for example the discussion in Ref. [21]), 
it follows that the SM fermions must be very weakly coupled to the new dynamics, with the 
exception of the top quark. 



The constraints on the dipole operators of Eq. (2.4) come from the current experimental 



limits on electric dipole moments (EDMs) and anomalous magnetic moments. The bounds 
on the neutron and mercury EDMs for example strongly constrain the dipole operators with 
u and d quarks. By using the formulas of Ref. 23 we find, with 95% probability, that: 



-7.01 x 10~ 6 < lm(c uB + c uW ) < 7.86 x 10" 
-9.42 x 10~ 7 < lm(c dB - c dW ) < 8.40 x 10" 



•1.62 x 10~ 6 < Im(c uG ) < 2.01 x 10" 6 , 



(2.17) 



-7.71 x 10~ 7 < Im(c dG ) < 5.70 x 10~ 7 , 
where the coefficients are evaluated at the low-energy scale /i ~ 1 GeV. According to the 



naive estimate (2.9), for 0(1) CP-violating phases these results imply a bound on {v/f) 2 at 
the level of 10~ 3 . In natural extensions of the SM, such a strong limit clearly points to the 
need of a symmetry protection mechanism. For a discussion, see for example Ref. 



21 



for 



the case of composite Higgs theories, and Ref. 24 for the case of SUSY theories 



Among the heavier quarks the most interesting bounds are those on dipole operators 



with top quarks 25 . These come from the experimental limit on the neutron EDM, 



- 1.39 x 10~ 4 < lm(c tG ) < 1.21 x 10" 4 



(2.18) 



the b — > 57 and b — > sl + l rates, 



0.057 < Re(ct W + c tB ) - 2.65 lm(c tw + c tB ) < 0.20 



and the tt cross sections measured at the Tevatron and LHC, 



(2.19) 



6.12 x 10~ 3 < Re(c tG ) < 1.94 x 10" 



(2.20) 
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All these bounds have 95% probability and have been derived by making use of the formulas 
reported in Ref. {25 1. Q It is worth noting that the bounds of Eqs. (2.19) and (2.20) are still 
about one order of magnitude weaker than the size of c t c (Hw an d c t B expected from the 
naive estimate (2.9) with (v/f) 2 ~ 0.1. Additional weaker constraints arise from the limits 



on anomalous top interactions based on top decays and single top production. From the 



results of Ref. 26 we find that, with 95% probability: 



1.2 < Re(c bW ) < 1.1 



-0.01 < Re(c tw ) < 0.02. 



(2.21) 



where the coefficients are evaluated at the scale /i ~ rut- 



in the lepton sector, the current measurements and SM predictions of the muon 27 28 



and electron 29 , 30 anomalous magnetic moments and the limits on their EDMs 31 , 32 



imply the following 95% probability bounds: 

-1.64 x 10~ 2 < Re{c eB - c eW ) < 3.37 x 10" 3 
1.88 x 10~ 4 < Re{c^ B ~ c^ w ) < 6.43 x 10~ 4 , 



(2.22) 



(2.23) 



-2.97 x 10~ 7 < Im(c eB - c eW ) < 4.51 x 10~ 7 , 

-0.26 < Im(c MB - c„ w ) < 0.29 , 

where the coefficients are evaluated at the relevant low-energy scale. Notice that the non- 
vanishing value of Re(c^B ~ c^w) follows from the known ~ 3.5<r anomaly in the (g — 2) of 
the muon (see Ref. [27] for an updated review). Among the bounds of Eqs. ( 2.21[ ), (2.22), 



(2.23) only those on Im(c e £ — c e w) and Re(c M £ — c^w) have the sensitivity to probe the 



values naively expected for these coefficients as reported in Eq. (2.9). In particular, the first 
one sets an upper bound on (v/f) 2 of order 10~ 3 for an 0(1) CP phase. 



3 Estimates of the effects on physics observables 



While the Lagrangian AC = ACsilh + ACf 1 + ACp 2 is completely general, the basis of 



operators of Eqs. (2.2)-(2.4) is particularly useful to characterize the interactions of the 



7 The coefficients are evaluated at the following scales: /i = m t (Eqs. (2.181 and (2.20)), /i = 

(Eq. wm). 
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Higgs sector. In fact, as already anticipated, one of the main results of Ref. |i] is that 
of identifying which operators, hence which observables, are sensitive to the strength of 
the Higgs interactions, rather than merely to the value of the New Physics scale M. In 
what follows we will discuss this point in greater detail and, starting from the analysis of 
Refs. (4 33 , we will try to highlight a possible strategy to determine whether the dynamics 



behind the electroweak breaking is weak or strong. Our analysis will be based on the naive 
estimates of the Wilson coefficients at the matching scale. In the next Section, we will discuss 
how the running from the matching scale to the weak scale affects these estimates. 

3.1 Operators sensitive to a strongly-interacting Higgs boson 

Let us start by considering the effects of the operators Oh, Ot, U) d,i and 0§: they modify 
the tree-level couplings of the Higgs boson to fermions, vector bosons and to itself. In the 
unitary gauge and upon canonical normalization of the Higgs kinetic term, the Lagrangian 



reads 34 



C = Id^hd^h - \m\h 2 - c 3 \ ) h* + . . . 



2 p 2 " 6 V v 

+ m 2 w W;W-»(l + 2c w ^ + ... I +-m z z Z^\ I + 2c z - + . . . J (3.24) 



h \ 1 2 _ _„ / . _ h 



I, — ,, r] ] V ' / 



where the Higgs couplings Ci=w,z,^,3, have been defined such that q = 1 in the SM, and 



v is defined by Eq. (2.5). Their expressions as functions of the coefficients of the effective 



Lagrangian (2.2) are given in Table [T] The shifts from the SM value are of order 



„2 2 2 

Hence, measuring the Higgs couplings probes the strength of its interactions to the new 
dynamics. Notice that the effective description given by AC neglects higher powers of (H/ f), 
and is thus valid only if the shifts in the Higgs couplings are small: 5ci ~ {v/f) 2 1. If the 
Higgs doublet is the NG boson of a spontaneously broken symmetry Q — > H, on the other 
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Higgs couplings 



ACsiLH 



MCHM4 MCHM5 



cw 



cz 



Cip (V> = u, d, I) 



c 3 



yy 



-77 



Cz~> 



cww 



czz 



CwdW 



CzdZ 



Czd-y 



1 - c H /2 
1 - c H /2 - 2c T 
1 - (c ff /2 + c^) 
1 + c 6 - 3c H /2 
8(« s /a 2 )c s 
8 sin 2 ^ c 7 



7 U' 









(chb — chw — 8 c 7 sm 2 8 w ) tan 6> 

-2 chw 
-2 (cfl"iy + c H b tan 2 9 w — 4c 7 tan 2 ^ sin 2 ^) 
-2(c w + chw) 
-2(c w + chw) - 2 (c B + c~hb) tan 2 6> w 

2 (cb + chb — c~w — chw) tan W 













-z 


Vi 


-z 


1 - 


2£ 


Vi 




1 - 


2^ 


Vi 





Table 1: The second column reports the values of the Higgs couplings Cj defined in Eq. (3.46) in 
terms of the coefficients Cj of the effective Lagrangian ACsilh- The last two columns show the 
predictions of the MCHM4 and MCHM5 models in terms of £ = (v/ f) 2 ; the effects of the heavy 
resonances have been neglected for simplicity, so that only the couplings cw,z,ip,3 are non- vanishing. 



The auxiliary parameter «2 is defined by Eq. (314B) 



hand, it is possible to resum all powers of (H/ f) by making use of the invariance under (non- 
linear) Q transformations. Such an improved effective Lagrangian thus relies only on the 
expansion in the number of derivatives. For example, in models based on the SO(h) / SO(&) 
coset [7 35 the couplings of the Higgs boson to W and Z are predicted to be cw = cz = 
cy = \/l — £, where £ = (v/f) 2 . The couplings to fermions, on the other hand, are not 
uniquely fixed by the choice of the coset, but depend on how the SM fermions are coupled to 
the strong dynamics. The last two columns of Table [l] report the predictions of the Minimal 
Composite Higgs Model MCHM4 [7] and MCHM5 [35], where the SM fermions couple linearly 
to composite operators transforming as the spinorial and fundamental representations of 
5*0(5), respectively. For simplicity, the predictions are derived by including only the effects 
of the Higgs non-linearities, and neglecting those from the heavy resonances, hence only 
the coefficients Cy, and C3 are non- vanishing. The models MCHM4 and MCHM5 will be 
considered as benchmarks in the rest of this work. 

In general, a shift of the tree-level Higgs couplings of order {v/f) 2 implies that the theory 
gets strongly coupled at energies ~ 47r/, unless new weakly-coupled physics states set in to 
regulate the energy growth of the scattering amplitudes. The dominant effect comes from 
the energy growth of the Vl^l ~* VlVl {V = W ± , Z°) scattering amplitudes, which become 
non-perturbative at the scale A s = Atcv/ ^/\ch\- A modified coupling to the top quark leads 
instead to strong VlVl —> ti scattering at energies of order A s = 1Qtc 2 v 2 / (m t \J\c u + c#|). 
The scale of New Physics is thus required to lie below, or at, such ultimate range of validity 
of the effective theory: M < A s . 



3.2 Operators sensitive to the scale of New Physics 

The operators Ow, Ob can be generated at tree-level by the exchange of heavy particles, for 
example heavy spin-1 states. In the unitary gauge they are written in terms of the following 
three operators [] 

(D»W+)W-"h, (d^Z^Z'h, {d^ v )Z v h (3.26) 

8 Here and in the following, derivatives acting on operators in the unitary gauge are covariant under local 
U(l) em transformations. Operators like {d^Z^^fh or (d^j^^ti obviously cannot be generated since 
they break the U(l) em local symmetry. 
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plus terms with zero or two Higgs fields. The fact that there are three possible operators 
in the unitary gauge indicates that their coefficients are related by one identity if the Higgs 



boson belongs to an SU(2) doublet, see Eq. (3.48). We will discuss this point in greater 
detail in Section 13.61 

It is easy to see that Ow, Ob give corrections to the tree-level Higgs couplings and 
generate quartic interactions with one vector boson and two SM fermions that contribute 
to the three-body decays h — > VV* — > Vipip. f\ Indeed, by making use of the equations of 
motion, Fj 

iD»W; u = gH^-^ u H - ig^- lv ^ , id^B^ = 9 - > - ig' jYj^ , (3.27) 
one can rewrite Ow arid Ob as 

O w = -2 H + ^ 2 {H^H)\D,H\ 2 + 0' Hq + 0' HL 
B = 2tan 2 ^ (-0 T + Y H *) , 



(3.28) 
(3.29) 



where the linear combination Ojj^, has been defined in Eq. (2.7). Upon the field redefinition 
H — > H — 2cw(H Jf H)H/v 2 , the operator [H^H)\D^H\ 2 can be rewritten in terms of those 
in Eq. (Q. Specifically, Eq. ( |3.28[ ) becomes: 



O 



w 



-60 H + 2 (O u + O d + O l )-%0 6 + 0' Hq + 0' HL . 



(3.31) 



From the estimates of c^y, cb and ch, ct, ty, c 6 in Eq. (2.9) one can see that the shifts to the 
tree-level Higgs couplings due to Ow, Ob are of order (mw/M) 2 , hence subdominant in the 
case of a strongly interacting Higgs boson. Notice that the couplings of the Higgs boson to 



9 We thank Riccardo Rattazzi for pointing this out to us 

10 For simplicity we have left a sum over all fermion representations ip understood in Eq. (3.271 
iii 



By means of Eqs. (3.291 and (3.31) it is thus always possible to remove Ow and Ob provided the 



coefficients of the other operators are shifted as follows: Cj — > c,; + Aq, with 

Ach = — 6cw, Act = — 2tan 2 9 'w c-b , Acq = ~8cw, Ac^—lcw 

a c' h = A? HL = c w 



(3.30) 



6 Ac Hq = - Ac Hu = -3 Ac H d = -2 Ac^l = -Ac m = -2 t'Av?0 w c B 



15 



W and Z get different shifts from Ob (since Act 7^ 0). In practice, the constraint (2.12) 
bounds this custodial-symmetry breaking effect down to an unobservable level, unless some 
fine tuning is in place in the combination c\y + cb so that cb can be large. Notice that 
despite the operator Or is generated after using the equations of motion, its contribution 



to Aei (corresponding to a non- vanishing T parameter 14 12 ) is exactly canceled by the 



vertex correction implied by the linear combination of fermionic operators which is also 



generated. _ This is of course expected, since Ow, Ob only contribute to €3, and not to e±. 

In general, the contribution of Ow, Ob to inclusive observables, in particular to the Higgs 
decay rates, is of order {m w /M 2 )\ 

5T(h ->■ VV) 

o w ,o B 



O(^). (3.32, 



T(h -> VV) 

where in this case VV = W^W*,Z^*'Z*,Z^*'j,jj. This implies that these operators are 
sensitive only to the value of the scale of New Physics M, and do not probe the coupling 



strength g*. From the quantitative side, the constraint (2.12) suggests that their effects in 
inclusive Higgs decay rates is too small to be observable. For example, we find that for small 
cw,b the tree-level correction to the WW and ZZ partial rates is well approximated by: p*| 

T(h -> W&W*) „ 00 Yih^Z^Z*) 1 nn/ 2n , 

nLwvw)L " 1 + 22 ° w ' r(/Lz(->z-)i " 1 + 2 -° {cw + tan ° w Cb) • 

(3.33) 

Notice that despite its custodial invariance, the operator Ow affects in a slightly different 
way the decay of the Higgs boson into WW and ZZ, due to the fact that at least one of the 
two final vector bosons is off-shell, m At the one-loop level Ow also contributes to the Higgs 



12 See for example Eq. (9.10) of Ref. 

13 Here and in the following our approximated formulas have been obtained by setting = 125 GeV and 
fully including the QCD corrections. The EW corrections to the decay rates are instead not included, since 
their effect on the numerical prefactor appearing in front of the coefficients q is of order (v 2 / f 2 )(a2/^) and 
thus beyond the accuracy of our computation. See Appendix 
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E.6 



for more details. 



It is easy to check that for > 2mz and on-shell decays one has: 

' ~l + 4c w , — 7^ ~ 1 + 4 gff + tan 2 g w c B ) . 3.34 

T(h -> WW) S m ro ZZ) SM K ' 

These formulas coincide with those of Eqs. (79)-(80) of Ref. j4j, which are thus valid only for on-shell decays. 
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decays into Z7 and 77 (while Ob does not). We find: 

T(h -> Z7) r(/i ^77) , 

r(A -> ^7)sm r(ft -> 77)sAf 

which agree with Eqs. (82) and (83) of Ref. [Ij. [^j For cv^b ~ 10 -3 the above approximate 
formulas imply corrections too small to be observed at the LHC. On the other hand, one 
could try to take advantage of the different predictions in terms of angular and invariant 
mass distributions which are implied by the dimension-6 operators compared to the tree-level 
SM prediction. The most promising strategy could be in fact that based on the analysis of 



the angular distributions of the final fermions 36 38 . In the ideal case in which one is able 
to kill completely the SM tree-level contribution by means of appropriate kinematic cuts, 
the relative effect of NP becomes of order 

dT(h^VV) 1 ( dY{h^VV) \ 16tt 2 

dn / { dn ) SM ~ 1 + CW < B > ( 3 - 36 ) 

which might leave room for observable effects even for cwb ~ 0(1CT 3 ). Clearly, a more 



precise assessment of the efficiency of such a strategy requires a dedicated analysis 39 



3.3 Operators generated at the one-loop level 

Let us now focus on the operators Ohw, Ohb, 7 and O g , which are generated at the one-loop 
level. In the unitary gauge, Ohw,hb,i are rewritten in terms of 

W+W-^h, Z^Z^h, w^h, Z^h (3.37) 

plus other terms with zero or two Higgs fields. Since the coefficients of the above four 



operators are functions of chw, c hb and c 7 , they are related by one identity, see Eq. (3.47). 



We will discuss this point in greater detail in Section 3.6 



As implied from the naive estimates (2.9), the contribution of Ohw,hb,i to the WW and 



ZZ inclusive rates is of order (VV = WW, ZZ) 



5T(h -> VV) 



T(h -> VV) 



0~,,Ohw ,Ohb 



rn 



2 

w 



167T 2 / 2 



15 



The easiest way to compute the one-loop contribution of Ow to the Z7 and 77 rates is by using Eq. ( 3.28 ) 



to rewrite this operator in terms of the others. Among the operators generated in this way, only Oh gives 



a contribution. Notice that if Eq. (3.31) is used instead, one has to take into account also the contribution 



of (O u + Od + Oi) and the shift to the Fermi constant induced by 0' Hq + 0' HL . 
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Although such an effect depends on the Higgs interaction strength, it is suppressed compared 



to Eq. (3.32) by a loop factor. We find that the following approximate formulas hold 
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r(h -). w(*)w*) SM 

T(h -> ZWZ*) 



1 + 3.7c#iy 



1 + 3.0 (chw + tan 2 #vy ch\b) — 0.26 c 7 



(3.40) 



r(/i -> z(*)z*) SM 

While the contribution due to chb and c 7 explicitly violates the custodial symmetry and 
thus differentiates WW from ZZ, the different numerical factor multiplying chw i n the two 
formulas above is due to the off-shellness of at least one of the two vector bosons, similarly to 



Eq. (3.33). Although there is currently no stringent bound on the coefficients chw,hb,~/, 

the 



estimate (2.9) suggests that their correction to inclusive rates is unobservable at the LHC. 
As discussed in the previous section, on the other hand, a study of the angular and invariant 
mass distributions of these decays can potentially uncover the effect of New Physics. In 



particular, an estimate similar to that of Eq. (3.36) can be derived also for chw,hb,j- 

The processes h —> 77, h — > Z7 and h — >■ gg (or equivalently gg —> h) can in principle 
test the Higgs interaction strengths much more powerfully, since they arise at the one-loop 
level in the SM. Naively one expects: 



r(/i -> ##,77,^7) 

16 For rrih > 2mz and on-shell decays, we find instead 
T(h -> WW) 



o 



Og,0~f,OHw ,Ohb 



T 2 



(3.41) 



T(h -> WW) SM 
T(h^ ZZ) 



1 



(3.39) 



1 + 8 (chw + tan 2 0^ c^s) — 16 tan 2 6w sin 2 9w c 7 . 



r(h -> zz) SM 

Comparing with the analog formulas in Eqs. (79) and (80) of Ref. [4j, we find that in these latter there is a 
missing factor 2 and the term proportional to c 7 was not included either. Notice also that the effect of the 



off-shellness of the gauge bosons is rather large, as one can see by comparing Eq. (3.39) with Eq. (3.40). 
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We find that the following approximate formulas hold to good accuracy for small q's: 

T(h ->■ qq) Att 
V ^ l + 22.2c 9 — 



T(h -)> gg) S M "2 

T(h ->■ 77) 47r . 

' 1 0.54 c 7 , (3.42) 



r(/i -» 77) 5 a/ a. 



cm 



1 + 0.19 (c w - c^b + 8c 7 sin 2 # w ) 



where we have conveniently defined 

_ V2G F m 2 w 

a 2 = — , (3.43) 

7T 

and by a em we indicate the value of the running electromagnetic coupling a em (q 2 = 0) in 
the Thomson limit. If the Higgs boson is a NG boson, the coefficients c g and c 7 are further 



suppressed by a factor (g@/g*) 2 , see Eq. (2.10), where g@ is a weak coupling. This implies 



that in this class of theories the corrections to T(h — » 77) and T(/i — )■ gg) depend only on 
the scale of New Physics and not on the Higgs interaction strength. In fact, in the case 
of minimal models with linear couplings, like for example the MCHM4 and MCHM5, the 



low energy theorem 40 , 41 implies that the leading contribution to the 77 and gg decay 



rates from the virtual exchange of heavy fermions is additionally suppressed 42 -45 due to a 



cancellation between the effect parametrized by c„ 7 and the one that follows from the shift 



in the top Yukawa coupling due to c u and ch (see Ref. 44 for an interesting exception). In 
general, in theories with a pNGB Higgs boson the local corrections to the rates T{h — > 77) 
and T(h — > gg) from 7 and O g are expected to be small and subdominant compared to the 
effect from the modified tree-level Higgs couplings. 

3.4 Fermionic operators 

The fermionic operators in A£p 1 are sensitive to the strength of the couplings of the Higgs 
boson and of the SM fermions to the new dynamics. They lead to contact corrections to the 
three-body decays h — > VV* — » Vipip which are naively of order 



5T(h -> Vjnjj) Q (^K 



Ci ■ (3-44) 
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Compared to the corrections from Ow and Ob, the effect of the fermionic operators is 
potentially enhanced by a factor (X^/g 2 ). In practice, the possibility of large fermionic 
couplings X^p is strongly constrained by LEP, see Eqs. (2.14)-(2.16). The only scenarios in 
which a large degree of compositeness of either the left- or right-handed quarks is not ruled 
out are those in which the corresponding operators in A£p 1 are not generated as due to some 



protecting symmetry (see for example Refs. 21,46,47]). Large corrections to the inclusive 
rate of the three-body decays h — > Vipip from ACpx are thus excluded, while the possibility 
of detecting the effects of these operators through the analysis of differential distributions 
should be explored, similarly to what has been discussed for Ow and Ob- 

Among the dipole operators in A£p 2 , those with light fermions are already strongly 
constrained by current precision data, but potentially sizable effects could still come from 
the operators involving the top quark. For example, the contribution of Ota to gg — > h, 
gg ti, gg — > tih is of order E 2 / (Wn 2 f 2 ) , where E is the energy scale relevant in the 
process. More in detail 



8v{gg -> h) 



CtG 



— 7 TTT" ~ CtG 

a(gg -> tt) m t 



$<j{gg — > tth) 



m 



o{gg -> h) 

:e we have ( 

that the experimental limit on the neutron EDM puts an upper bound on the imaginary 



where we have defined c tG = Re(c tG ) (mj/m^) ~ m 2 t / \16tt 2 f 2 ) ~ 3 x 10 3 (v 2 / f 2 ). Notice 



a(gg tth) 

r 2{2\ 



CtG 



2 • 



(3.45) 



part of CtG at the 10 level, see Eq. (2.18), which indicates that this is currently the most 
sensitive experiment on Im(c t c)- Some mechanism is however required to suppress the imag- 
inary parts of the dipole operators involving light fermions, in order to satisfy the stringent 



constraints of Eq. (2.17). By the same mechanism also Im(ctcO could be suppressed, so that 



the processes of Eq. (3.45) are essential to probe the contribution of O t c due to Re(ctc). 



From Eq. (3.45) and the naive estimate of c t c it follows that the most sensitive process is 
perhaps gg — > ti, in particular the events at large invariant mass, although a precision larger 
than the one currently achieved is required to constrain (v/f). To this aim, the analysis of 



differential distributions and spin correlations could be a successful strategy 48 25 . The 
NP contribution to the process gg — > tih can in principle get the largest enhancement from 



a cut on y/s, but the small rate might limit the actual sensitivity achievable at the LHC 49 



Finally, additional information comes from the experimental limits on top anomalous cou- 
plings obtained at the Tevatron and the LHC, although their sensitivity on NP is expected 
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to be much smaller by naive estimate. The operator O t w, m particular, gives the largest 
effect and generates the anomalous coupling gR{g/mw)bL^ v W~ v tR 26 . Naively one ex- 
pects gn = (Am t /m w ) c t w ~ rnt m w / i^® 712 f 2 ) = 1.5 x 1CT 3 (v/ f) 2 , an effect too small to be 
observed even for / of order v. 



3.5 Non-linear Lagrangian for a Higgs-like scalar 

Summarizing, by working in the unitary gauge and in the basis of fermion mass eigenstates, 
the effective Lagrangian relevant for Higgs physics reads as follows |34| 



£ = \dji d"h - -mlh 2 - c 3 - 
2 M 2 h 6 



3m 2 



h 6 - 



— „ rl 1 \ 



ip=u,d,l 



+ m 2 w W+W 



h 

l + 2c w - + .. 

v 



h 



+ -m 2 z Z„Z» 1 + 2c z - + . . . + 



(3.46) 



+ [c ww W+W~^ 



czz 



^77 



liv _|_ 99 ga ^yafiu 



+ (c WdW (w-D^w^ + h.c.) + c zdz z u d,z^ + c zdl z u d^ u ) \ + ■ ■ ■ 



where, we recall, v is defined in Eq. (2.5). We have shown only terms involving up to 



three bosonic fields, and we have omitted in particular those involving fermions that follow 



from A£p x + ACp 2 . Their form can be easily derived from Eqs. (2.3) and (2.4). The 



relations between the couplings appearing in Eq. ( |3.46 ) and the coefficients of the dimension- 
6 operators in Eq. (2.2) are reported in TableJT} It is worth noting that the same Lagrangian 



(3.46) applies also to the case in which the electroweak symmetry SU(2)l x U(1)y is non- 
linearly realized and h is a generic CP-even scalar, singlet of the custodial symmetry, not 



necessarily connected with the EW symmetry breaking. Indeed, each of the terms in (3.46) 



being invariant under local U(l) em transformations, can be dressed up with the Nambu- 
Goldstone bosons that are eaten to form the longitudinal W and Z polarizations and made 
manifestly SU(2)l x U(1)y gauge invariant [50] (see also Ref. [51] ). The explicit expression 



in such a basis has been given in Refs. 52,53 at the level of four-derivative operators. In this 



sense the effective Lagrangian (3.46) is a generic tool to understand the origin of the newly 



discovered boson and the role it plays in the electroweak symmetry breaking dynamics. It 



is valid for arbitrary values of the couplings Cj appearing in Eq. (3.46), and it can be used 
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to make computations of observable quantities at a given order in an expansion in E/M 
and in asul^-, where by the latter we indicate the generic SM loop expansion parameter. 
That is in full analogy with other well-known effective theories, see Ref. [3l. It should be 



stressed that, according to a well established methodology and similarly to Eq. (2.2), in 



this effective Lagrangian all quantum fluctuations associated to short-length modes (high- 
energy modes) have already been considered and are parametrized by local operators with 
an increasing number of derivatives, while quantum fluctuations (loop diagrams) involving 
the light modes still have to be taken into account. For instance, top loops will give an 



additional contribution to the on-shell /i-gluon-gluon coupling. While Eq. (3.46) is general 



the effective Lagrangian (2.2) assumes that h is part of an SU(2)l doublet and further relies 
on the expansion in powers of H/f. As such, it is valid only in the limit of small deviations 
of the Higgs couplings from their SM values and up to corrections of order 0(v 2 /f 2 ). 



3.6 Implications of custodial symmetry 



Another difference between the non-linear Lagrangian (3.46) and the SILH Lagrangian (2.2) 



is that the first one contains two more free parameters. This means that there are two 



relations among the couplings of Eq. (3.46) which hold at the level of dimension-6 operators 



if the Higgs is part of a doublet. As noticed in Sections (3.2) and (3.3), the first identity 



relates cww, czz, c^ 7 and c 77 , while the second relates cwaw, czaz and cza-y- They read: 



cww — czz cos 2 9w = cz~{ sin 29 w + c 77 sin 2 ^ 

Czd-y 



cwaw — CzdZ cos 9 W 



sin 29w ■ 



(3.47) 
(3.48) 



In fact both identities are not special to the case in which the Higgs is a doublet, but are a 
general consequence of custodial symmetry. This latter is accidental in the SILH Lagrangian 
if one restricts to the operators that lead to derivative couplings of the Higgs to vector 
bosons. Starting at the dimension-8 order, it is possible to write cutodial-breaking operators 



that lead to couplings that violate the relations (3.47) and (3.48). For instance 

,2 



csww 9 
m 2 w v 2 



(H^W° u a a H) (H^W b » u a b H) + (H^a a H) {D^W^f (h^Ph) (3.49) 
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gives rise to 



czaz = — 4c 8 h/, cza-y = —4: tan 9wC8\v , 

Czz = 8 COS^OwCswW 5 c Zy = 4 sin 29\yCgww j C-yy = &SlIl 2 9wC8WW 



(3.50) 



and the relations (3.47) and (3.48) are not fulfilled. 

A third relation holds on the non-derivative couplings cw and cz if one assumes that 



custodial symmetry is an invariance of the Lagrangian (2.2), so that ct = 0; it reads: 



c w = c z . 



(3.5i; 



As said above, while all three identities (3.47), (3.48) and (3.51) are a consequence of custo- 



dial symmetry, the first two are accidental at the level of dimension-6 operators if the Higgs 
is part of a doublet. 



To show that Eqs. (3.47), (3.48) and (3.51) follow from custodial invariance, let us con- 



sider the case in which the EWSB dynamics has a global SU(2)l x SU(2) r symmetry, and 
imagine to fully gauge this group by enlarging the hypercharge to a whole triplet of SU (2) R . 
In this case the diagonal custodial SU(2)v is exact even though g' ^ g. The left and right 
gauge fields couple to the conserved currents of SU (2) L x SU (2) R and the interactions among 
two gauge fields and the Higgs boson are fully characterized in momentum space by three 
form factors: 



Fll)%(pi,P2) L^Ujh. + (T LR )Z( Pl ,p 2 ) L\B? v h + (r RR )%( Pl ,p 2 ) KJVJi . 



'v 



(3.52) 



Here pi, p 2 are the momenta of the gauge fields and each form factor can be computed in 
terms of a Green function with two conserved currents, = (JfJ^\h). In addition to 
the usual massive W and Z bosons, which form a triplet of the custodial group, in this 
case there is a whole triplet of massless SU(2)v gauge fields (the photon plus its charged 
companion), VZ. The mass eigenstates and are related to the left and right gauge fields 
through a rotation by an angle 8\y, where tan#y^ = g'/g. Their cubic interactions with the 
Higgs boson are thus characterized by three form factors, which are linear combinations of 

U = 



The two operators in (3.49) give rise to the oblique parameter U, see for instance Ref. 



12 



-c S w - Zcshw while S = c 8 hw- 
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those in Eq. (3.52): 
T V v 

1 vv 

F „ . 
1 VV 



sm 2 9 w Y LL + 



sin2# 



(Y LR + Y RL ) + cos 2 6»vk Tij/j 



sin20w , 2 . 2 , sm26 w 
— 2 — r LL + (cos 6T LjR - sin fcT^J — T RR 

cos 6V T LL (T LR + T RL ) + sin 6>n/ r RjR , 



(3.53) 



where we have defined r^(pi,p 2 ) = r^(p 2 ,Pi)- Notice, in particular, that in this case the 
same form factor T vv describes the interaction of two W's and two Z's to the Higgs boson, 
as due to custodial invariance. 

The physical limit where only SU(2)l x U(1)y is gauged is obtained by simply switching 
off the unphysical R 1 ' 2 fields. The interactions of two neutral vector bosons to the Higgs are 



still described by the relations of Eq. (3.53), where T zz 



F- - F 

1 VVi 77 



VV 



and T 



VV 



In the charged sector, instead, the W corresponds to a pure left gauge field, since it has no 
mixing with right-handed ones. This implies that its form factor is given by the last formula 



of Eq. (3.53) with 6w = 0, that is: Tww = Tll- The four physical form factors are linear 



combinations of the three defined in Eq. (3.52), and are thus related by one identity: 

iw(pi,P2) -r^(pi,p 2 )cos 2 ^ = (r^(pi,p 2 ) + r^(p 2 ,Pi)) sm ^ w + T^(p u p 2 ) sm 2 e w . 

(3.54) 

Notice that this relation is a consequence of our initial assumption of SU(2)l x SU(2) r 
invariance of the EWSB dynamics. The custodial SU(2)y is broken in this case only by 



the gauging of hypercharge. For g' = the custodial symmetry is unbroken and Eq. (3.54) 
implies T ww = T zz . It is straightforward to derive the relations (3.47), (3.48) and (3.51) 



from Eq. (3.54). At quadratic order in the momenta, the form factors can be computed from 
the effective Lagrangian ( |3.46[ ); one has: 



r£(pi,P2) 



2m 2 w c w rT ~ 2c ww P% - c wdw {P? v + P 2 ^) 
2m 2 z c z rf v - 2c zz P% - c zdz (P? u + If) 

— 2C Z ~ / Pi2 — CZQjP^ 



(3.55) 



- 2r P^ u 



where we have defined P^ = rj^p 2 — PxP\, P^ = ?7 M!/ p 2 — P2P2 an d P\2 = V^Pi'Pz ~ P1P2 
This is in fact the most general decomposition which follows at the 0(p 2 ) level for an on-shell 
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Higgs boson by assuming CP invariance and requiring that: i) the IViy, Yzz and T 77 form 
factors are symmetric under the exchange {pi, /i} -H- {p 2 , v}; ii) the T 77 and Y Zl form factors 
satisfy the Ward identities implied by U(l) em local invariance: 

Pi^( Pl ,p 2 ) = = p 2v Y%( Pl ,p 2 ) , p 2v Y%(pu P2 ) = . (3.56) 

Additional structures proportional to pi^ and p 2v can be omitted since they give vanishing 
contributions both when the vector bosons are on-shell and when they decay into a pair of 



fermions by coupling to the corresponding conserved current. Inserting Eq. (3.55) into (3.54) 



one then obtains the identities (3.47), (3.48) and (3.51) 



From the above discussion it follows that if custodial symmetry is an invariance of the 



EWSB dynamics, the effective Lagrangians (3.46) and (2.2) have the same number of free 



parameters, in terms of which all observables can be computed. This is true also if one 



considers the fermionic operators (for a Higgs doublet these are listed in Eqs. (2.3) and 



(2.4)), as long as one focuses on terms with one Higgs boson. This means that by using 
single-Higgs processes alone, one cannot distinguish the case in which the Higgs boson is 
part of a doublet from the more general situation. The only possible strategy to this aim 
is exploiting the connection among processes with zero, one and two Higgs bosons which 



is implied by the Lagrangian (2.1) at 0(v / f ) and does not hold in the case of the more 



general non-linear Lagrangian. As a consequence of such connection, the bounds that EW 
and flavor data set on operators with zero Higgs fields severely constrain the size of the NP 



effects in Higgs processes, as discussed in Section (2.1). If one were to find that single-Higgs 



processes violate these constraints, this would be an indication that the Higgs is not part of 
a doublet. Furthermore, processes with double Higgs boson production can be predicted to 
a certain extent in terms of single-Higgs couplings, and can thus be used to probe the nature 



of the Higgs boson 54 



4 Implementing the Higgs effective Lagrangian beyond 
the tree level 



In this section we address a few issues related to the use of the effective Lagrangians (2.1) 



and (3.46) beyond the tree level, as required to make Higgs precision physics without as- 
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suming the validity of the Standard Model. While the methodology is well established and 
various examples of its application exist in several different contexts, we think that a dedi- 



cated discussion can be useful to better clarify some specific points (see also Ref. 55 for a 
recent discussion). As an illustrative though important example, we will consider the cal- 
culation of the Higgs partial decay widths, and show how the corrections from dimension-6 
operators can be incorporated in a consistent way. As a by-product of our analysis and to 
better demonstrate its applicability, we present in Appendix IE] a modified version of the pro- 



gram HDECAY 56 which features a full implementation of the effective Lagrangian ACsilh, 



Eq. (2.2), as well as its generalization to the case of a non-linearly realized EW symmetry, 



Eq. (3.46). 



A first difficulty which arises when using either Eq. (2.1) or (3.46) is the presence of mul- 
tiple expansion parameters. For generic values of the Higgs couplings Cj, the validity of the 
effective Lagrangian ( 3.46[ ) is based on a double perturbative expansion in the SM couplings, 



q;5m/47t, and in powers of E/M. The effective Lagrangian (2.1) further assumes (v/f) <C 1, 
which implies small shifts in the Higgs couplings: q = l + ^Q, with 5ci < 0(v 2 /f 2 ). All these 
expansion parameters must be properly taken into account when performing calculations. 
Furthermore, the non-renormalizability of the effective theory implies the presence of addi- 
tional divergences compared to the SM case which must be absorbed by a renormalization 
of the Wilson coefficients of local operators. 

4.1 RG evolution of the Wilson coefficients 

Let us discuss the issue of the renormalization and RG evolution of the Wilson coefficients 
first. As done in the previous sections, we will assume that the Higgs boson is part of an 



SU(2)l doublet and use the Lagrangian (2.1). Since we are only interested in the divergent 



structure of the diagrams, it is convenient to work in the limit of unbroken SU(2)l x U(1)y 
and compute the Green functions in terms of the Higgs doublet H. The only 1-loop dia- 
grams which generate additional logarithmic divergences are those featuring one insertion 
of the effective vertices from dimension-6 operators. By dimensional analysis, further inser- 
tions of the effective vertices lead to power-divergent contributions to dimension-6 opera- 
tors (which are irrelevant to determine the RG running) and log-divergent contributions to 
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higher- dimensional operators. The same counting holds also at higher loop level: the only 
log-divergent contribution to dimension-6 operators comes from diagrams with one inser- 
tion of the effective couplings, and is thus suppressed by extra powers of the SM expansion 
parameter asul^- This is in analogy with the renormalization of the pion effective La- 
grangian in the chiral limit, see Ref. [57]. It thus follows that the RG equation is linear 
and homogeneous in the Cj, and different operators with the same quantum numbers will in 
general mix with each other. At leading order in asM, with asM = &em, «2,a s , respectively, 
in the case of electromagnetic, weak and QCD corrections, one has 

*G0 = (% + 7<°> ^ k*(£)) c,(M) , (4.57) 

where 7^ is the leading-order coefficient of the anomalous dimension. Some elements of the 



„(o) 
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anomalous dimension matrix 7^ have been recently computed in Refs. 

In the case in which the Higgs boson and possibly the SM quarks (in particular the top 
and the bottom) are strongly coupled to the new dynamics, the leading RG running effect 
comes from loops of these particles and can be as large as Ac$/ci(M) ~ (g 2 /16n 2 ) log(M//i) 
or (A 2 /16tt 2 ) log(M///). This must be compared to the effects of order (g| M /167T 2 ) log(M//i) 
from loops of gauge fields. For example, the insertion of c B in the diagram (a) of Fig. [I] leads 
to a renormalization of Ow+b = Ow + Ob'- 

c w+b (ii) = c w+B (M) - ^ g log (^) c H (M) , (4.58) 

where «2 has been defined in Eq. (3.43). It is well known that this RG running is associated 
with the IR contribution to the e 3 parameter, and the same coefficient 7w+b h = — 1/6 



can indeed be extracted from self-energy diagrams 60 . From the estimates of Eq. (2.9), 
c H (M) ~ 0(v 2 /f 2 ), c WjB (M) ~ Oimly/M 2 ), it follows that the correction to c w+B from 
its RG evolution down to the scale \i is of order Acw+b/c\v+b(M) ~ (g 2 /16ii 2 ) log(M//x) 
as anticipated. Similarly, the insertion of ch^ into a loop of fermions, like in diagram (b) of 
Fig. [TJ leads to a renormalization of cw and c B : 

Ac W)B « N c ^ log (^) c H ^M) , (4.59) 

where = 3 is a color factor. In this case the RG correction is of order (A^/lGvr 2 ) log(M//i) 
compared to the UV value of the coefficients, as one can immediately verify by using the 



estimates (2.9). 
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Figure 1: One-loop diagrams relevant for the RG running of cw and cb- Dashed, continuous 
and wiggly lines denote, respectively, a weak doublet H, a fermion and a vector field V = W, B. 
The symbol ® denotes the insertion of the effective vertex from Oh (in diagram (a)) or Oh*/; (in 
diagram (b)). 




(a) (b) 



Figure 2: One-loop diagrams relevant for the RG running of ct- Dashed, continuous and wiggly 
lines denote, respectively, a weak doublet H, a fermion and a hypercharge field B. The symbol ® 
denotes the insertion of the effective vertex from Oh (in diagram (a)) or On^> (in diagram (b)). 



Loops of EW gauge fields give corrections which are suppressed by a weak loop factor 
(g 2 /167r 2 ), and the associated RG evolution is therefore generically small. An important 
exception is the case in which the Wilson coefficient has a value suppressed at the scale M. 
For example, if the dynamics behind the EW symmetry breaking is custodially invariant, then 
ct(M) = 0. The insertion of ch into a loop of hypercharge gauge bosons, as in diagram (a) 
of Fig. [2j renormalizes ct and gives 

5r(M) = \ tan 2 ** g log (±) c H (M) . (4.60) 

Compared to the naive estimate of Eq. (2.9), ct{M) ~ 0(v 2 /f 2 ), valid in absence of custo- 
dial symmetry, the above correction is further suppressed by a factor (g /2 /167r 2 ) log(M//z). 
Although small, such a low-energy value of ct has a strong impact on the EW precision tests 
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performed at LEP |60|. 18 On the other hand, it is too small to be observable through a 
measurement of the Higgs couplings at the LHC A similar renormalization of ct also follows 
from loops of SM fermions through the insertion of cjj^p, as illustrated by diagram (b) of 
Fig. [2] The explicit calculation for the case of a composite right- and left-handed top quark 
was performed for example in Ref. [61] . Naively, the effect goes like 

AdT ~ N ^ l0g (M)~ CH ^ M) > (461) 
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and is of order (y^/g'YiX^/g*) 2 compared to the one from loops of hypercharge. 

In general, although small, the RG evolution of the Wilson coefficients due to EW loops 
must be properly taken into account in order to precisely match the experimental results 
obtained at low energy with the theory predictions at high energy. This is even more true 
in the case of QCD loop corrections, which can be large and will affect the coefficients of 



the dimension-6 operators with quarks and gluon fields. The effect of the running of 



the Wilson coefficients can be easily incorporated in programs for the automatic calculation 



of production cross sections and decay rates by using the effective Lagrangian (2.1) and 
identifying the coefficients appearing there as their values at the relevant low-energy scale. 
This is, in particular, how we proceeded to build eHDECAY, the modified version of HDECAY 
presented in Appendix |E} 

4.2 Decay rates at the loop level with the effective Lagrangian 

In addition to the short-distance effects discussed above, which are parametrized in terms 
of the evolution of the coefficients of local operators, one-loop diagrams also lead to long- 



18 For example, c T (m z ) ~ 1CT 3 for c H (M) ~ 0.1. 

19 Notice that in case of a sizable fermion coupling A^, a numerically larger contribution to ct comes from 

fermionic loops with two insertions of Ch^. The corresponding diagram is quadratically divergent, so that 

it gives a threshold correction to c"t at the scale M, but does not contribute to its running. An explicit 

calculation can be found in Ref. [6l] for the case of a composite top quark. Naively the effect is of order 

Act ~ N c (v / f) 2 (\^/16tt 2 )(X^ /g*) 2 , and can be numerically large. For example, if both and couple 

with the same strength X tL = X tR ~ \Jg*Vt to the new dynamics, then it follows Act ~ N c (v/ f) 2 (y 2 /IQtt 2 ). 
20 Notice that g 2 c g is not renormalized at one-loop by QCD corrections. This follows from the RG- 

invariance of the operator (f3(g s ) / g s )G M „G Ml ' which contributes to the trace of the energy-momentum ten- 



sor 62 64 . See also the recent discussion in Ref. 158 . 
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distance corrections to the observables under consideration. Specifically, while short-distance 
effects are related to the divergent terms, the long-distance contributions correspond to 
the finite parts and are defined in a given renormalization scheme. In general, the decay 
amplitude can be expanded as follows: 



A = Aq + A i + AA + AAi + ... (4.62) 

where Aq M (Af AI ) is the tree-level (one-loop) SM amplitude, and AA (AAi) is the tree- 
level (one-loop) contribution from the dimension-6 operators of the effective Lagrangian in 



Eqs. (2.2)-(2.4). The dots denote higher-loop contributions as well as the corrections due to 
higher-order operators. 

Let us consider for example the decay h — > W^*'W*. In this case the operators that can 
contribute at tree-level are Oh, Ow, Ohw, O^w, 0' H ^, as well as Ouud in the case in which 



the off-shell W decays into a pair of quarks. Based on the naive estimates of Eq. (2.9) and 
according to the discussion of Section [3j we can quantify the various effects encoded by AAq 
as follows: 

^fv 2 A u A d \ ^ (v 2 \\\ A / Em^ 

+ c H ud x O — — — +c m xO\—^\ + c^w x O f 



(4.63) 

Here E = is the relevant energy of the process and we have conveniently defined each of 



21 In the strict sense this equation is valid for the genuine EW corrections only, while for simplicity we 
include the (IR-divergent) virtual QED corrections to the SM amplitude in the same way. The corresponding 
real photon radiation contributions to the decay rates are treated in terms of a linear novel contribution 
to the Higgs coupling for the squared amplitude in order to obtain an infrared finite result. Pure QED 
corrections factorize as QCD corrections in general so that their amplitudes scale with the modified Higgs 
couplings. However, they cannot be separated from the genuine EW corrections in a simple way. 
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the 0(1) parameters q to be equal to Ci{rrih) divided by its naive estimate in Eq. (2.9): 

. M 2 , 4 



-^Ci{m h ), 
v 



H,TM, 



m 



-Ci{m h ) 



W,B 



ir 



16n 2 f 2 _ 



rn 



-Ci{m h 
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C; 



gU 2 
\ 2 



2 A(m h ), 



gl P 



(4.64) 



CHud 



CHud(m h ) 



A U A<2 u 2 

When the Higgs boson is pNGB, the two parameters c g and c 7 are not of order one but are 
further suppressed by a factor g^jg 2 . From Eq. (4.63) one can see that the contribution of 
the dipole operators O^w is suppressed by {m^/m h ) compared to that of Ohwi while that 
of Ojjud an d 0' H ^ is expected to be small given the existing constraints on the couplings 



A^, (see the discussion in Section 2.1). The dominant NP contribution thus comes from the 



terms in the first line of Eq. (4.63), among which the one proportional to ch is the leading 



effect for > g. The 1-loop electroweak amplitude Af M gives a contribution of order 
Af M /A^ M ~ (a 2 /47r). We thus see explicitly that AA and Af M encode the NLO corrections 
in the three expansion parameters which we are considering: a2/47r (electroweak expansion), 
E 2 /M 2 (derivative expansion) and v 2 /f 2 . The contribution due to 1-loop diagrams with one 
insertion of the effective vertices has not been computed yet, but we can easily estimate its 



size: 



AAi 



+ 



(4.65) 



where the dots denote the subleading terms due to the other operators. The terms shown 



in Eq. (4.65) arise from the same 1-loop diagrams that give the SM amplitude Af M , where 
each of the Higgs couplings gets shifted by c#, c u and c 6 . By neglecting the unknown AAi 
one is omitting terms of order (w 2 // 2 )(a 2 /47r), that is, of the same size of the tree-level 
contribution due to the operator Ohw, see Eq. (4.63), since E = rrih ~ mw- This latter 



contribution can be easily computed and it is included in the formula of the decay rate to 
WW (and similarly that of Ohw an d Ohb to ZZ is also included) implemented by the 



program eHDECAY discussed in the Appendix (see Eq.(E.131) and also Eqs.(3.40), (E.138) 



(E.139)). The addition of the tree-level correction from Ohw is clearly the first step towards 
a full inclusion of the 0[(w 2 // 2 )(a2/47r)] corrections, where the missing part will have to be 
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computed from 1-loop diagrams featuring one insertion of Oh-, O u and 0§. It is worth noting 
that these diagrams, in general, contain logarithmic divergences which must be reabsorbed 
by a renormalization of the Wilson coefficients and contribute to their RG evolution as 
explained in the previous section. The finite part is the contribution to AAi which awaits 
to be computed. 

By approximating the amplitude as A ~ Aq + Af M + AA one obtains the following 
formula for the decay rate: 2 ^_ 

2 



T(W^W*) = T S M (W^W*){ 1 + Re[(A s Q M Y (Al M + AA )] 

[ \ A I 

+ 0! 




(4.66) 



where Tq M (W^W*) denotes the tree-level SM decay rate. For simplicity, we have not 
shown terms involving powers of E 2 /M 2 among the neglected contributions, since for E = 
rrih ~ m\y one has E 2 /M 2 < v 2 /f 2 if g* > g. As mentioned, this formula incorporates 
the 0(v 2 /f 2 ), 0(ct2/47r) and 0(m 2 h /M 2 ) corrections (NLO in the perturbative expansion), 
and can be easily implemented in existing codes for the automatic computation of the decay 
rate. The inclusion of the 0(m 2 h /M 2 ) tree-level correction due to Ow is justified as long as 
g* < Att, since it is parametrically larger than the neglected 0[(v 2 /f 2 ) (0:2/47?)] terms by a 
factor (167r 2 /g^). Notice that in the limit of large deviations of the Higgs couplings from their 
SM values, (v/f) 2 ~ 0(1), the neglected terms of 0[(v 2 / f 2 )(a 2 /ATr)] become as important 
as those included through Af . In order words, a proper inclusion of the EW corrections in 
the limit v ~ / requires a complete 1-loop calculation where each of the diagrams is rescaled 
by the appropriate coupling factor. 

A similar discussion applies to the Higgs decay into a pair of fermions, h — > x/ji/j. In this 
case only Oh and (ip = u, d, I) contribute at tree level, 

£p<W)=(^ + a,)xo(£), (4.67) 
while the one-loop EW diagrams featuring one effective vertex give a correction of order 

^ m ^ H xo(^) + c,xo(^) + «xo(^) + ... (4.68) 



A$ MKYYJ n \P Air) ' w \P 4tt J 1 "° " ~ \f2 4?r 



2 The same remark as in footnote 
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applies. 
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where the dots indicate the subleading terms due to the other operators. The calculation 
of AAi has not been performed yet, while the 1-loop EW corrections are known in the SM, 
Af M . Their inclusion is thus possible as long as (v/f) <C 1, so that the neglected terms 
in AAi are subleading. The case of QCD radiative corrections is different, since at leading 
order they factorize with respect to the expansion in the number of derivative and fields and 
can thus be resummed up to higher orders. In the case of the Higgs decay into a pair of 
quarks one can for example approximate A ~ Aq M + Af + AAq and obtain the following 
formula for the decay rate: 



2:i 



T(qq) = r s Q M (qq) k qcd j 1 + J a Jm^ Re [ (^o M ) * «< + AA )} 



+ 




(4.69) 



where ^Q M {qq) is the SM tree-level rate and k® cd encodes the QCD corrections. This for- 
mula includes the leading 0(v 2 /f 2 ), 0(a2/4vr) and QCD corrections. Mixed electroweak 
and QCD corrections have been also included by assuming that they factorize, as the non- 



factorizable terms are known to be small. Compared to the decay rate into WW, Eq. (4.69) 
apparently does not include corrections of order m|/M 2 . While there is indeed no operator 
whose contribution starts at that order, such corrections can arise from subleading contri- 
butions to ch and c^. For example, the tree-level exchange of heavy fermions can lead to a 
wave-function renormalization of the SM ones, which can be re-expressed in our notation as 
a contribution to of order \ 2 h v 2 /M 2 . 

A similar resummation of the QCD corrections also works for the decay h — )■ gg. In this 
case the SM tree-level amplitude vanishes, A$ M = 0, while the leading contribution arises 
from the 1-loop exchange of top quarks. The two-loop EW corrections are known in the SM 
and give a correction of order Af M /Af M ~ a2/47r. Among the dimension-6 operators, only 
O g contributes at tree-level, 

AAo (gg) = c 9 xo( v ^] . (4.70) 



A SMW 9 \f 



As discussed in Section^ (see Eq. (2.10)), the above estimate is suppressed by an additional 



3 The same remark as in footnote 
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applies. 
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factor (g%/g 2 ) in the case of a NG Higgs boson. At the one-loop level one has 



AAl to) = ( t + x o ( t )+ a» x o ( v l 4-, I • (^i) 



Thus, the one-loop effect of 0# and O u is expected to be as important as the tree-level 



one from O g , and even larger if the Higgs is a NG boson, as discussed in Section 3.3 This 
is in fact not surprising, since c g arises at the 1-loop level in minimally coupled theories, 
while ch and c u can be generated at tree level. The contribution from the dipole operator 
O t G is suppressed by a factor y 2 /16ii 2 compared to that from Oh and O u , as expected from 
the fact that etc is generated at the 1-loop level in minimally coupled theories. For this 
reason it will be neglected in the following. It should be noted that without a complete 
computation of the NLO EW corrections of order (a2/47r)(t; 2 // 2 ), the LHC data on Higgs 



physics are not sensitive to the range of values of c t c expected using the naive estimate (2.9) 
with {v/f) 2 ~ 0.1. Furthermore, we stress that in order to distinguish the effect of Ota from 
that of O g , the tth channel should be measured (49] (single top production in association 
with the Higgs could also provide complementary information [65]). Also in this case, there 
are no operators giving m 2 h /M 2 corrections, although these terms will in general appear as 
subleading contributions to c g , ch and c u , as discussed above. It is well known that higher- 
order a s corrections are large, so they must be included consistently in our perturbative 
expansion. This can be done easily in the approximation rrih <C 2m t , which is reasonably 
accurate for rrih = 125 GeV. In such a limit one can integrate out the top quark and match its 
one-loop contribution to that of the local operator O g . Then it trivially follows that the QCD 
corrections associated to the virtual exchange and real emissions of gluons and light quarks 
below the scale m t factorize in the rate, the multiplicative factor being the same for both 
the top quark and New Physics terms. By approximating A ~ Af M + Af A/ + AAq + AA\, 
one arrives at the following formula for the h — )• gg decay rate: 

r(gg) = rf M ( gg ) n soft {c 2 eff + M(Af M Y ( A l Mc *ff + AA ° + AA ^ c e //)] 

(4.72) 

/ / r~ \ / rl.i \ / M .-> \ - \ 1 





An f 2 

where Tf M (gg) is the 1-loop SM decay width. The factor c e // includes all the depen- 
dence on m t and accounts for virtual QCD corrections to Af M above that scale, while 
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K so ft parametrizes the soft radiative effects. By using Eq. (4.72), the existing four- loop 
66 69 and K so f t 70 -72 allow one to include the QCD corrections up 



calculations of c e ff 
to N 3 LO. 

The contributions to the decay h — > 77 follow a similar pattern as for h — > gg. At tree 
level: 

AA 



A SM 



(77) 



£ 7 XC \ p 



(4.73) 



At one loop: 



(77) =c H xO 



P 

+ c^iy x O 



+ c u x O 



P 



+ c w x O 



rn 



n? 



+ (c tw + c tB ) x O 



M 2 
2 



V 



(4.74) 



16tt 2 / 2 / V./ 2 l ^ 2 , 

The 2-loop electroweak corrections have been computed in the SM and can be included for 
(y 2 /f 2 ) 1, so that unknown 0[(f 2 // 2 )(a2/47r)] effects arising from 2-loop diagrams with 



one effective vertex are negligible. From Eq. (4.74) one can see that the 1-loop contribution 



due to Ohw is of the same order as such neglected terms. Furthermore, it is currently not 
known and thus not included in the formulas implemented in eHDECAY, see Appendix [Ej The 
1-loop correction from Ow, on the contrary, is parametrically larger by a factor (167r 2 /g 2 ) 
and should be included for g* < 4tt. The easiest way to compute it is by rewriting Ow 
in terms of the other operators through the equations of motions |4], see Eq. (3.31). The 
result is finite and has been included in the formulas implemented in eHDECAY. The 1-loop 
correction due to the dipole operators is suppressed by a factor y 2 /16n 2 and can be neglected. 
Approximating A ~ Af M + v4f M + AA + AAi one finds: 



T(77) = rf'(77) |l + |^ Re[(Af M )* (A™ + AA + AAj] 
+ 0\ 




(4.75) 
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AA °(Z 1 )=c,xo(^)+(c HW -c HB )xo(^) , (4.76) 



Finally, the estimate of the corrections to h — > •yZ is the following: 
a sm v" u ~ -7 -~ i ^ ) + ( Cw ~ CHfi ) x lj2 

^7) = «-xo(J) +4xog) +a^xo(^) + iwX o(^ 

' v ( 4 - 77 ) 




P 16tt 2 y V/ 2 167r2 . 

In this case the 1-loop electroweak corrections are not known in the SM, so that the formula 
for the decay rate reads: 

T(Z 7 ) = rf»(Z 7 ) {i+i^RbIK")* (AA + AA,)] +o(( J)' , 

where only the contributions from 0#, O u and Ow should be retained in AA\ for consistency. 

Through the above discussion we sketched how the effective Lagrangian can be imple- 
mented beyond the tree level in the calculation of physical quantities. In the case of the 
Higgs partial decay widths, in particular, we have seen how the EW and QCD corrections 
can be included consistently with the expansion in the number of fields and derivatives. As a 
more concrete illustration of these considerations, we have written a modified version of the 
program HDECAY, which we dub eHDECAY, where the corrections from all the local operators 



of the effective Lagrangian (2.2) are included at NLO. A detailed description of the code is 
given in Appendix [Ej where more explicit formulas for each of the Higgs partial widths are 
provided. 



5 Discussion 

The discovery of a resonance with a mass around 125 GeV similar to the long-sought Standard 
Model Higgs boson brings the exploration of the electroweak symmetry breaking sector under 
quantitative scrutiny. The LHC experiments, together with those at the Tevatron, report 
the signal strengths, i.e. the product of the Higgs production cross section times its decay 
branching ratio, for various final state channels. The main task of the community is now to 
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interpret these data and understand the implications for the theory of New Physics that is 
expected to lie beyond the weak scale. 

The EW oblique parameters provide a bound on the scale of New Physics but do not 
give detailed information about the nature of the NP sector. In order to understand how the 
weak scale is stabilized at the quantum level, i.e. how the hierarchy problem is solved, one 
crucial question is whether EW symmetry breaking proceeds by weak or strong dynamics. 
The direct observation of new degrees of freedom would provide a straightforward answer. 
But a glimpse of New Physics can also be caught from a dedicated study of the Higgs boson 
itself, and in particular from a measurement of its couplings, if a departure from the SM 
predictions is ever observed. It is useful to parametrize the deviations from the SM by the 



effective Lagrangian of Eq. (2.1). By measuring its Wilson coefficients q one can infer what 
kind of UV theory completes the SM. 

If the coupling strength of the Higgs boson to the NP sector is of the order of the SM weak 



couplings, <?* ~ g, then our power counting (2.9) shows that the coefficients of the operators 
that can be generated at tree-level, Oh,O u ^,i,Ow and Ob, are expected to be all of the 
same order, m^/M 2 , where M is the typical mass scale of the NP spectrum, unless some 
special selection rule suppresses some of them. It is instructive to examine the predictions 
of the archetypal example of weakly-coupled UV completions: the Minimal Supersymmetric 
Standard Model (MSSM). First, i?-parity protects the EW oblique parameters from any 
tree-level contributions, hence cw and cb are of order (m^/M 2 )^/^^) and thus small. 
Second, the couplings of the lightest Higgs boson to the massive gauge bosons are given by 
cy = sin(/3 — a), where a is the rotation angle to diagonalize the CP-even mass matrix and 
tan /3 is the ratio of the vacuum expectation values of the two neutral CP-even Higgs bosons. 
In the decoupling limit, a — > (3 — n/2, one has cy = l + 0{m 4L z /m\ { ) 1 where m# is the mass of 
the heaviest CP-even scalar (for a general treatment of the decoupling limit see for example 
Ref. |73|). This means that at tree-level the deviations of the Higgs-gauge boson couplings 



are generated by dimension-8 operators 74 , while ch arises only through loop effects and is 
naively of order (m^/M 2 ) (a 2 /47r). At the same time, the couplings to up- and down-type 



37 



quarks read, respectively, 

Cu = + = 1 + 2 cos^/3 cos 2/3 + 



cos a 


= 1 + 2 


m| 


sin /3 




sin a 


= 1-2 


m 2 z 


cos /3 


m H 



4 

m H 



c d = = 1 - 2 — f sin^/3 cos 2/3 + 



(5.79) 



For moderately large tan/3 this implies q ~ m 2 z /m 2 H , while c u is further suppressed by a 



factor ~ 1/ tan 2 /3 (see for example Refs. 75 76 and the recent discussion in Ref. [T7j ) . A 



pattern with small values of ch, cw, cb and c u but with a ~ 15% enhancement of the Higgs 
coupling to down-type quarks due to q, for example, would be indicative of the MSSM with 
large tan /3 and the additional Higgs bosons around 300 GeV. Generic two-Higgs doublet 
models lead to a similar pattern of couplings, while models where the Higgs boson mixes 
with a scalar that is singlet under the SM gauge group can generate ch at the tree level. In the 
MSSM, loops of light stops or staus as well as charginos can also give sizable contributions to 
the effective couplings of the light Higgs boson to photons and gluons, with c g , c 7 satisfying 



the naive estimates (2.9). For example, loops of stops lead to c g ~ (g 2 /167i 2 )(m 2 v /m 2 ), 
where = y t or A t /mi. 

This situation has to be contrasted with the case of strongly coupled theories. There, our 



power counting (2.9) singles out cr, c U) d as the dominant Wilson coefficients (cq controls only 
the Higgs self-interaction and measuring it at the LHC will be challenging), while cw and cb 
are suppressed by the ratio (g/g*) 2 - Furthermore, a composite Higgs boson can be naturally 
light if it is the pseudo Nambu-Goldstone boson associated to the dynamical breaking of a 
global symmetry of the strong dynamics. This implies that the coefficients c g and c 7 will also 
be suppressed by a factor (g$/g*) 2 , where g@ is some weak spurion breaking the Goldstone 
symmetry. The modifications in the gluon-fusion production cross section and in the decay 
rate to photons are thus controlled by ch and c u . 

The harvest of data collected by the LHC certainly calls for a definite theoretical frame- 
work to describe the Higgs-like resonance and compute production and decay rates accurately 
in perturbation theory without restricting to the SM hypothesis. Effective Lagrangians are 
one of the tools at our disposal to achieve this goal. Elaborating on the operator classifi- 
cation of Ref. jl), we estimated the present bounds on the Wilson coefficients and provided 
accurate expressions for the Higgs decay rates including various effects that were previously 
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omitted in the literature. Assuming that the observed Higgs-like resonance is a spin-0 and 
CP-even particle, we discussed two general formulations of the effective Lagrangian, one of 
which relies on the linear realization of SU(2)l X U(l)y at high energies. One of the ques- 
tions that can be addressed by considering these two parametrizations is whether the theory 
of New Physics flows to the SM in the infrared, that is, whether the Higgs-like resonance 
is part of an EW doublet. If all the Higgs signal strengths measured at the LHC converge 
towards the SM prediction, it would be a very suggestive indication that indeed the Higgs 
boson combines together with the longitudinal components of the W and Z to form an EW 
doublet, since any other alternative requires some tuning to fake the SM rates. On the other 
hand, the doublet nature of the Higgs boson would be less obvious to establish if the signal 
strengths exhibit deviations from their SM predictions. We have pointed out that, if the 
EWSB dynamics is custodially symmetric, it is not possible to test whether the Higgs boson 
is part of a doublet by means of single-Higgs processes alone. A direct proof can come only 
from processes with multi-Higgs bosons in the final states [54] , which are however challenging 
to study at the LHC. Precisely establishing the CP nature of the Higgs boson is another 
question that also requires accurate computations. If there is little doubt that the observed 
resonance has a large CP-even component, the possibility of a small mixing with a CP-odd 
component remains alive, and dedicated analyses will have to be performed to bound the 
mixing angle between the two components. To this aim too, an effective Lagrangian includ- 
ing the CP-odd operators listed in Appendix [C] provides the theoretical framework where 
this question can be addressed quantitatively. 

The absence so far of direct signals of New Physics at the LHC indicates that the road 
to unveil the origin of the electroweak symmetry breaking might be long and go through 
precision analyses rather than copious production of new particles. For such a task, the well 
established technology of effective field theories is the most powerful and general tool we have 
to analyze the Higgs data and put them into a coherent picture together with the existing 
experimental information without assuming the validity of the Standard Model. There is 
still time for the LHC to disprove this pessimistic eventuality by reporting the discovery of 
new light particles or large shifts in some of the Higgs couplings. It is clear, however, that if 
the New Physics continues to remain elusive, a precise investigation of the Higgs properties 
will become the most urgent programme in high-energy physics both for the experimental 
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and the theoretical community. 
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A SM Lagrangian: notations and conventions 

In this Appendix, we collect the conventions used throughout this paper. The field content 
decomposes under SU(3)c x SU{2) L x U(l)y as 

H = (1, 2, 1/2), L\ = (1,2, -1/2), r R = (1, 1, -1), (A.80) 
qi = (3,2,1/6), i4= (3,1,2/3), <f R = (3, 1, -1/3), (A.81) 

where the hypercharge is defined as Y = Q — T^l, and i — 1, 2, 3 is a flavor index. The action 
of the gauge group is fully characterized by the conventions used to define the covariant 
derivative. For instance, for the left-handed quark doublet, we have 

D,q L = (d, - Z -g s X a g; - Z -ga l W; - % -g'B^\ q L (A.82) 

where A a , a = 1 ... 8, and a 1 , i = 1 ... 3, are the usual Gell-Mann and Pauli matrices. Ac- 
cordingly, the gauge-field strengths are defined as 

G% = d»9u - d,g; + g s f abc 9l9t , (A.83) 
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where f abc are the SU(3) structure constants. 

The Yukawa interactions of the up-type quarks involve the Higgs charge-conjugate dou- 
blet defined as 

H c = ia 2 H*. (A.84) 
The renormalizable Lagrangian of the SM thus reads: 

C SM = - \G%G a ^ - \w; v W^ - \b, v B^ + {D»H)\D^H) 

+ i {L L YD^L L + IrYD^u + qirfDyflL + u R YD,u R + djffDpdR) ( A - 85 ) 
+ /4jH^H - \(H ] Hf + (y u q L H c u R + y d q L Hd R + y t L L Hl R + h.c.) 

B Electroweak Chiral Lagrangian in non-unitary gauge 

We report here the expression of the EW chiral Lagrangian valid in a generic gauge and in 
the most general case in which the SU(2)i x U(l)y is non-linearly realized. For simplicity, 
we will restrict to the case in which the EWSB dynamics has a custodial invariance. The 
scalar h is assumed to be CP-even and a singlet of the custodial symmetry, and does not 
necessarily belong to an SU(2)l doublet. The Lagrangian can be expanded in terms with 
an increasing number of derivatives 

C = C + C EW sb , Cewsb = -V(h) + £ (2) + £ (4) + . . . (B.86) 

where £ contains the kinetic terms of the SU (3) c x SU (2) L x U (l)y gauge fields and of the 
SM fermions, Cewsb describes the sector responsible for EWSB, and V(h) is the potential 
for h [34]: 

V(h) = + c 3 I (M) e + . . . (B.87) 

Under the request of SU(2)y custodial symmetry, the longitudinal W and Z polarizations 
correspond to the NG bosons of the global coset SU (2) L x SU {2) R / SU (2)y and are described 
by the 2x2 matrix 

E(x) = exp (ia a X a (x)/v) , (B.88) 
where a a are the Pauli matrices. SU(2)l x U(1)y (local) transformations read as 

£(x) -> U L E(x)U Y , U L = exp(ia a L a a ) , U Y = exp(ia Y a 3 ) (B.89) 
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and the covariant derivative is defined by 



= c^E - % -j-W; <x a E + % -^B, E a 3 



(B.90) 



At the level of two derivatives one has (34|: 



£ (2) = \{d,hf + ^Tr (D^LTE) (l + 2c v ^ + 



V2 

v 

V 



(0 



A§ (ug ) ,dg ) )S(0,dg J )""(l + ^^ + 



+ h.c. 
+ h.c. 



(B.91) 



-^(*f)s(o,i«) T (i + ^ + ••-)+/, 



c. 



where the dots stand for terms with two or more Higgs fields and an implicit sum over flavor 
indices i,j = 1, 2, 3 has been understood. After rotating to the fermion mass eigenbasis and 



by choosing the unitary gauge E(x) = 1, the sum of (B.87) and (B.91) coincides with the 



first two lines of Eq. (3.46) with cw = cz = cy. 



At the level of four derivatives, there are 6 independent bosonic operators 24 which affect 
cubic vertices with one h field: F3 



= c' ww W; v W^ a h - + c' WB Tr (St W« v o a E B^a 3 ) £ + c' BB B^ h - 



+ ^- D»W a au Tr (tfa'ffivE) h - d^B. v Trf E+z^E 



m w 



a 3 U 



(B.93) 



2 - v 

The dots stand for terms which have two or more h fields or do not lead to cubic vertices, 
for the complete list of bosonic operators in & A >. In the unitary gauge, 



see Refs. 



52 



53 



24 The operator 2 = 0/™w) Tr[(-D/*£)t(A,E)] ^^/i that appeared in Eq. (B.85) of the first version 
of this paper is actually redundant and can be eliminated in terms of the operators 0' ww and 0' WB . For 



C2 and c'. 



wb 



-wb 



+ tan 2 #iy C2 can be used to remove C2- 



example, the shift c' H 

25 Another convenient basis, which can be more easily compared to Eq. (3.46), is one in which the first two 



operators of Eq. (B.93) are replaced by 



W a Tr 

v Y [lis 11 



Sfh, 



Bpi, Tr 



d v h. 



(B.92) 



This is in fact the basis adopted in Rcf. 
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Eq. (B.93) coincides with the last three lines of Eq. (3.46). More specifically, the coefficients 



cwwi c zzi cz 7 , c 77 can be written as linear combinations of c' ww , c' BB , c' WB , 

2d u 



'7? ^77 

c ww — * '-■WW 

c Z z = 2(cos 2 9 w c' ww - 2 sin 9 W cos 9 W c WB + sin 2 ^ c BB ) 

.2/ 



c 77 = 2(sm 2 9 w c' ww + 2 sin 6 W cos 6 W c' WB + cos 2 6 w c' BB ) 
c Zl = 2 (sin W cos 9w c' ww + cos c' WB — sin 9w cos 9w c' BB ) 
while cwawi c zaz can be expressed in terms of c' w , c' B : 

cwaw = 4:c' w 
czaz = 4c' w + 4 tan 9 W c' B 
cza-y = 4 tan 9 W c' w -Ad B . 



(B.94) 



(B.95) 



Notice that Eqs. (B.94) and (B.95) are directly implied by Eq. (3.53), which follows from 



custodial invariance. It is simple to verify that the identities (3.47) and (3.48) are satisfied 



by the couplings appearing on the left-hand sides of respectively Eq. (B.94) and (B.95). 



The above discussion shows explicitly that every operator in Eq. (3.46) can be dressed 
up with NG bosons and made manifestly invariant under local SU(2)l x U(1)y transforma- 
tions. P3 



The part of Eq. (B.86) which does not depend on the Higgs field h coincides with the 



non-linear chiral Lagrangian for SU(2)l x U(l)y [78], in the limit of exact custodial sym- 
metry. This latter assumption can be relaxed by specifying the sources of explicit breaking 
of the custodial symmetry, i. e. its spurions, in terms of which one can construct additional 
operators formally invariant under SU(2)l x U(1)y local transformations. For example, the 
list of operators that follows in the case in which custodial invariance is broken by a field 



with the EW quantum numbers of hypercharge has been recently discussed in Ref. 53 



Since the choice of quantum numbers of the spurions is model-dependent (and in fact the 

strongest effects are expected to arise from the breaking due to the top quark, rather than 

hypercharge), we do not report here any particular list of operators, and prefer to refer to 

the existing literature for further details. 

26 Notice that h is invariant under SU(2)l x SU(2)fi (hence SU(2)l x U(1)y) transformations. In the 
case in which h belongs to an SU(2)l doublet H, this follows from the fact that h parametrizes the norm of 
the doublet: WH = (v + h) 2 /2. 
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C Relaxing the CP-even hypothesis 



If one relaxes the hypothesis that h is CP-even, there are six extra dimension-6 operators 



that need to be added to the effective Lagrangian (2.2): 
ic H w 9 



AC 



CP 



rn 



(D^Hya l (D u H)W l 



ic H Bg 



w 



fJ,U 

2 



772 



(D^Hy(D u H)B i 



[IV 



IV 



+ H^HB^B^ + ^ H^HG a av G a ^ 
m w m w 

~ 3 ~ 3 

+ e ijk Wi V Wl P W k / + f abc G a» G ^p£ 



(C.96) 



IV 



c fj. 

p 



where the dual field strengths are defined as F^ = \e^p a FP a for F = W, B, G (e is the 
totally antisymmetric tensor normalized to e i23 = 1)- Furthermore, the coefficients of the 
operators involving fermions will be in general complex numbers. 

In the case of the effective chiral Lagrangian with SU(2)l x U(1)y non-linearly realized, 



there are four additional operators, to be added to those of Eq. (B.93), which can affect 
cubic vertices with one h field: 



AC 



(4) 
CP 



~c' ww W^W^ a 



h 



+ ~c' 



WB 



Tr 



h 

v 



h 



+ c' BB B, u B^-+ } 



c gg no. nauv 



(C.97) 



In the unitary gauge, both Lagrangians ACcp and AC^ P are matched onto: 



AC 



(4) 
CP 



c. 



99 no 



J fit/* 



\ 2 

+ (cwdw (W~D,W + ^ + h.c. 



Gt,G apu 



h 

v 



c zaz Z v d,Z^ 



5 ZdJ z„d fl r u )^ + --- 



(C.98) 



When the EW symmetry is linearly realized, the coefficients of Eq. (C.98) are related to 



the Wilson coefficients of Eq. (C.96) through the same relations reported in Table [T] with 
the simple exchange q — > Ci (and with cw = cb = 0). In the non-linear case, cwaw = 
czaz = czd~j — while cww, czz, c 77 and c^ 7 are given in terms of the Wilson coefficients of 



Eq. (C.97) by relations identical to the ones of Eq. (B.94) (with q — > Ci and C2 = 0). 

Finally, it should also be noted that when the CP-invariance assumption in the Higgs 
sector is relaxed, the couplings c u ^i are allowed to take some complex values. 
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D Current bounds on dimension-6 operators 



In this Appendix we explain how we derived the bounds on the coefficients of the dimension-6 



operators reported in Section 2T For a given observable we construct a likelihood for the 
coefficients c* as follows: 

L{ci) oc exp [-(Osm + SOfa) - O exp f/{2 A0 2 exp )} , (D.99) 

where O exp ± AO exp is the experimental value of the observable, Osm denotes its SM pre- 
diction and 50 (q) is the correction due to the effective operators. If several observables 
constrain the same coefficients Cj, the global likelihood is constructed by multiplying those 
of each observable. We include the theoretical uncertainty on the SM prediction by integrat- 
ing over a nuisance parameter whose distribution is appropriately chosen. We then quote 
the bound on a given coefficient by marginalizing over the remaining ones. 



Let us consider for example the bounds of Eqs. (2.14) and (2.15). To derive them we used 



the EW fit performed in Ref. 20 by the GFitter collaboration, and constructed a likelihood 
for the various coefficients by computing their contributions to the Z-pole observables. For 
the latter, we used the SM predictions and experimental inputs reported in Table 1 of 



Ref. 20 , treating the uncertainties on the SM predictions as normally distributed. We 
performed two separate fits: one on the coefficients of the operators involving the light 
quarks (u, d, s), and one on those with charged leptons and heavy quarks (c, b). We thus 
neglected, for simplicity, the correlations between these two sets of coefficients. The relevant 
observables in the first fit are T tot , a had and R[. They depend on the Wilson coefficients only 
through the following linear combination: 

I = (k ~ ^ sin 2 ^ (c Hq i + c' Hql ) + f-i + i sin 2 #^ (c Hql - c' Hql + c Hg2 - c' Ha2 ) 



+ \ sin 2 6»iy c Hu - - sin 2 9 w (c H d + c Hs ) 
6 o 



(D.100) 



which with 95% probability is constrained to lie in the interval 

- 0.62 x 10~ 3 <l< 1.2 x 10~ 3 . (D.101) 

Although there are no further observables at the Z-pole which can resolve the degener- 
acy implied by this result, we thought it useful to report the limits that one obtains from 
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Eq. (D.101) by turning on one coefficient at the time. These are the bounds reported in 



Eq. (2.14). 



The second fit, performed on the coefficients of the operators with leptons and heavy 
quarks, makes use of all the observables at the Z pole and counts 7 unknowns, specifically: 
Xi = {{c Hq 2 + c' Hq2 ), c Hc , (c Hq 3 - & Hq3 ), cm, cm, {chl + c' HL ), (c H l - c' HL )}. For simplicity we 
assume lepton universality, and thus take the coefficients chu c~hl, c' hl to be the same for all 
the leptonic generations. In terms of the above variables, the result of the fit is summarized 
by their central values Xi, standard deviations <Ji and by the correlation matrix py-: 

1-3 



- . -/ 

C-Hql T C Hq2 
CHc 



CHq3 - c' Hq3 



CHb 

chi 

chl + c'hl 
chl - c'hl 



(-5.5 ±5.1) x 10" 
(4.3 ± 11) x 1(T 3 
(3.0 ±3.0) x 10~ 3 
(-3.5 ± 1.5) x 10~ 2 
(9.6 ±8.6) x 10~ 4 
(-6.4 ± 15) x 10~ 4 
(-1.3 ±0.75) x 10~ 3 



(D.102) 
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The limits of Eq. (2.15) have been obtained by making use of the above formulas and 



marginalizing over all the coefficients except the one on which the bound is reported. 



For the limits of Eqs. (2.11) and (2.12) we have used the fit on S and T performed in 



Ref. 20 , by marginalizing on one parameter to extract the bound on the other. 



To derive Eq. (2.17) we have used the theoretical predictions of the EDM of the neutron 



and mercury given in Ref. 23 in terms of the dipole moments of the quarks (see Eqs. (2.12), 
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(3.65) and (3.71) of Ref. |23|), and the experimental results for these observables given 



respectively in Ref. 79 and Ref. 80 . We included the theoretical errors by assuming that 



they are uniformly distributed within the stated intervals. Only two linear combinations of 
the coefficients q can be constrained in this way, since two are the observables at disposal: 

[lm(c dG ) + 1.3lm(c dB - c dW )] - [Im(c nG ) - 0.64Im(c uB + c uW )} 

(D.104) 

''in . ''in , 

h 



h 



m 



2m n 2m d 
—j- Im(c nG ) H 2~ lm i C dG) ■ 



m 



m 



Using m u = 2.3 MeV and m d = 4.8 MeV we obtain, with 95% probability: 

-1.59 x KT^GeV- 1 <l x < 1.78 x 1(T 12 GeV^ 1 

-1.82 x 10" 12 GeV" 1 < l 2 < 1.37 x 10" 12 GeV" 1 . 
From the above result, by turning on one coefficient at the time, one obtains the limits 



(D.105) 



given in Eq. (2.17). The bound on lm(c t c) of Eq. (2.18) has been similarly derived from 



the neutron and mercury EDMs by following Ref. [25] and making use of the formulas given 
there. 



The limits of Eq. (2.22) have been obtained from the experimental measurements of the 



electron [29 j and muon (27) anomalous magnetic moments and their SM predictions (taken 



respectively from Ref. [30] and Refs. 27 28]). In this case we have included the theoretical 
errors by assuming that they are normally distributed. All the remaining bounds reported in 



Section 2.1, namely those of Eqs. (2.19)-(2.21 ) and Eq. (2.23) have been obtained by simply 



translating into our notation the results given in the references quoted in the text. 



E eHDECAY: a full implementation of the Higgs ef- 
fective Lagrangian into HDECAY 



The program HDECAY 56 was originally written for the automatic computation of the Higgs 



partial decay widths and branching ratios in the SM and in its Minimal Supersymmetric 
extension (MSSM). Its latest version 5.10 includes the possibility of specifying modified cou- 
plings for up-type quarks, down-type quarks, leptons and vector bosons in the parametriza- 



tion of Eq. (3.46), as well as of including the effective couplings c„„, c 77 and c^ 7 . We present 
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here a modified version of the program labelled eHDECAY. It is available at the following 



URL: http://www-itp.particle.uni-karlsruhe.de/~maggie/eHDECAY/, In addition to 
the features already present in HDECAY, the new program includes the effective couplings 



cww, czz, cwaw and czaz, and thus fully implements the non-linear Lagrangian (3.46). 
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In fact, similarly to HDECAY 5.10, it also includes the possibility of choosing different cou- 
plings of the Higgs boson to each of the up and down quark flavors and lepton flavors. In 
this sense the program does not assume neither custodial symmetry nor flavor alignment. 



As explained in the text, Eq. (3.46) describes a generic CP-even scalar h at NLO in the 
derivative expansion. If h forms an SU(2)i doublet together with the longitudinal polariza- 



tions of the W and the Z, the Lagrangian can be expanded as in Eq. (2.2) for (v/f) <C 1; 
in this case the values of the Higgs couplings c, are given in the second column of Table [TJ 
The program eHDECAY provides an option in its input file where the user can switch from 



the non-linear parametrization of Eq. (3.46) to that of the SILH Lagrangian Eq. (2.2). The 
user can also choose to set the values of the Higgs couplings to those predicted at tree-level 
in the benchmark composite Higgs models MCHM4 and MCHM5, see the last two columns 
of Table E 

Similarly to the original version of HDECAY, all the relevant QCD corrections are included. 
They generally factorize with respect to the expansion in the number of fields and derivatives 
of the effective Lagrangian, and can thus be straightforwardly included by making use of the 
existing SM computations. The inclusion of the EW corrections is less straightforward and 



can be currently done in a consistent way only in the framework of the Lagrangian (2.2) 
and up to higher orders in {v/f). Going beyond such approximations would require fully 
dedicated computations which at the moment are not available in the literature. In eHDECAY 
the user has the option to include the one-loop EW corrections to a given decay rate only 



if the parametrization of Eq. (2.2) has been chosen. The same EW scheme used by HDECAY, 
with Gf-i itlw and rriz taken as input parameters, is also adopted in eHDECAY. The sine of 
the Weinberg angle is defined as 



2 

HI 

sin 2 9 w = 1 - -f , (E.106) 



2 ' 

m z 



27 Notice that the operator proportional to cza-y does not affect the decay h — > Z7 as long as the photon 
is on-shell. 
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following the so-called on-shell scheme 81 . Derived quantities in this scheme are also the 



electromagnetic coupling and the weak coupling. To describe this latter we have conveniently 



defined the parameter a 2 in Eq. (3.43). The formulas implemented in the program are thus 
written in terms of only the input parameters or their derived quantities sin 6w and a<i- 
The only exception to this rule is given by the decay rates Y{h — > 77) and T(h — > Zj), 
where we use the experimental value of the electromagnetic coupling in the Thomson limit, 
oiemifj 1 = 0), in order to resum large logarithms for on-shell photons. 

Below it follows a detailed discussion of how the New Physics corrections are incorporated 
for each of the Higgs decay modes. We report explicitly the formulas implemented in the code 
and their level of approximation in the perturbative expansion of the effective Lagrangian. 
In all the following expressions, as explained in the text, the coefficients of the dimension-6 



operators of the SILH Lagrangian (2.2) and those of the derivative operators of Eq. (3.46) 



must be identified with their values at the relevant low-energy scale = m^. 



E.l Decays into quarks and leptons 



Upon adopting the effective description of the non-linear Lagrangian (3.46) and working at 
leading order in the derivative expansion, the Higgs boson partial decay width into a pair of 
fermions is obtained by rescaling the tree-level SM value rQ M (ipip) by a factor c 2 . The QCD 
corrections to the decay widths into quarks which are currently available for the SM case 



include fully massive NLO corrections near threshold 82 and massless 0(a>i) corrections far 



above threshold 83 -85 . Also, large logarithms can be resummed through the running of the 
quark masses and of the strong coupling constant. All these QCD corrections factorize with 
respect to the tree-level amplitude and can therefore be incorporated as done in HDECAY for 
the SM case. The decay rate can be written as follows: 

T(^)\ NL = 4lf^V>) [1 + 5i,K QCD ] , (E.107) 
where Tq M denotes the leading-order decay width, 8^ = 1(0) for ip = quark (lepton) and 



k qcd enco( j es th e QCD corrections. 28 This is the formula implemented by eHDECAY in the 

28 There is one caveat, however. In the case of decays into strange, charm or bottom quarks there are two- 
loop diagrams which involve loops of top quarks coupling to the Higgs boson. They need a rescaling different 
from c.0. It has been correctly taken into account by the appropriate modification factor c^Ct (ip — c : s : 
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case of the non-linear Lagrangian (3.46). It is valid up to corrections of 0(m 2 h / M 2 ) in the 
derivative expansion and of 0(a2/47r) from EW loops. These latter corrections are available 
in the SM but contrary to the QCD ones do not factorize. Their inclusion in the case 
of generic Higgs couplings thus requires a dedicated calculation, which is not available at 
present. 

In case of the SILH parametrization, where the deviations of the Higgs couplings from 
their SM values are assumed to be of 0(v 2 /f 2 ) and small, the decay rate can be written as 

2 



SILH 



r n 



c H - 2cw, + 



\A™\ 



Dp (a*SMaSM\ 



[1 + 5^k qcd ] , (E.108) 



where A^ M , Af^ w are, respectively, the tree-level and EW one-loop 86 amplitudes of the SM. 
In this case the one-loop EW corrections can be easily included if one neglects terms of 
0[(a2/4:7r)(v/ f) 2 ] (the same remark as in footnote 21 applies). In particular, mixed contri- 



butions up to 0[(o;2/47r)(a s /47r) 4 ] have been included by assuming that the electroweak and 
QCD corrections factorize, as the non-factorizable contributions are small. From the view- 



point of the expansion in inverse powers of the NP scale, formula (E.108) includes corrections 



of order 0(v 2 /f 2 ). It neglects terms of 0(v 4 /f 4 ) } 0[(a 2 /4vr)(t;//) 2 ], 0[(a 2 /4vr) 2 ]. 
E.2 Decay into gluons 



Upon selecting the Lagrangian (3.46), the rate into two gluons is computed in eHDECAY by 
means of the following formula: 

2 



GFOt 2 m\ 
[nl= 4^3 



q=t,b,c 



C "ff Ksoft 



+ 2Re( £ f^ /2 (V 27rC - 



\q=t,b,c 





27TC gg 


1 °eff Ksoft + 





Ksoft (E.109) 



+ 0" Yl C q A* 1/2 (T q )c q/ A 1/2 (T q/ )^ NLO (T q ,T q/ ) 



q,q'=t,t 



where r q = \mrjmr h and the loop function is defined as 

A1/2 (r) = \r [1 + (1 - r) / (r)] 



(E.110) 
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with 



f(r) 



arcsm 



In 



1 + y/l 



-, 2 



The first term corresponds to the one-loop contribution from the top, bottom and charm 
quarks, whose couplings to the Higgs boson are modified with respect to their SM values. 
The second and third terms encode the effect of the derivative interaction between the Higgs 
boson and two gluons generated by New Physics. Naively c gg ~ (a s /4:7r)(g 2 v 2 /M 2 ), so that 
the correction from the effective interaction can be as important as the one from the top 
quark if (g 2 v 2 /M 2 ) » 1. No expansion is thus possible in c gg in the general case. 
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r > 1 



r < 1. 



(E.lll) 



The QCD corrections have been included up to N 3 LO in Eq. (E.109) in the limit of 



heavy loop particle masses, similarly to what is done in HDECAY for the SM. In this limit the 
effect of soft radiation factorizes and is encoded by the coefficient K S oft- The coefficient c e ff, 
instead, takes into account the correction from the exchange of hard gluons and quarks with 
virtuality q 2 3> m 2 . More in detail, for rrih <C 2mj, one can integrate out the top quark and 



obtain the following five-flavour effective Lagrangian 

C cS = -2^Gfc, G° aflu G°rh 



(E.112) 



where bare fields are labeled by the superscript 0. The renormalized coefficient function C\ 
encodes the dependence on the top quark mass m t . The coefficients K so f t and c e // are thus 
defined as 

71 



K 



soft 



c eff 



lmU GG (q 2 = m 2 h ) 



2ml 



12ttCi 
ah\m h ) 



(E.113) 



where Il GG {q 2 ) is the vacuum polarization induced by the gluon operator. The N 3 LO ex- 
pression of the coefficient function C\ 66 69 in the on-shell scheme and that of Im I1 GG can 
be found in Ref. 72 . At NLO the expressions for n S oft and c e // take the well-known form 



,NLO 
v soft 



1 + 



a 



NLO 



71 



73 7 
T ~ 6' 



^NLO 
-eft 



1 + 



a 



NLO 



11 



71 



4 ' 



(E.114) 
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where here a^ LO is evaluated at the scale and computed at NLO for Np = 5 active 
flavours, while it has been consistently computed up to N 3 LO for the inclusion in eHDECAY. 



The last line in Eq. (E.109) contains the additional mass effects at NLO QCD [71] in the 
top and bottom loops, encoded in K (r g , T q >), which have been explicitly implemented in 
HDECAY and taken over in eHDECAY. While the mass effects for the top quark loops play only 
a minor role, below the percent level, for the bottom loop contribution the mass effects for a 
125 GeV Higgs boson amount to about 8% relative to the approximate NLO result. Hence, 



formula (E.109) includes the QCD corrections at N 3 LO (i.e. at 0(a h s ) in the decay rate), 
and neglects terms of 0{m\/ M 2 ) in the derivative expansion. 



When the SILH Lagrangian (2.2) is selected, on the other hand, eHDECAY computes the 



decay rate into gluons by means of the following approximate formula: 



SILH 



Gpoi 2 s m\ 
4\/27r 3 



\ Yl ( 1 - c H-c q -c q/ )A* 1/2 (T q/ )A 1/2 {T q )c 2 eff K soft 



9 

q,q'=t,b,c 



+ 2Re £ -A{ 



/2 V T <?J 



\q=t,b,c 



167T Cg 

a 2 



Ceff Ksoft 



q=t,b,c 



c eff ^soft 



+ 



9 (± - ch - Cq - Cq')A* 1/2 (T q ) A 1/2 (r q ,) K NL °(T q ,T q ,) 



q,q'=t,b 



(E.115) 

The last line contains the mass effects at NLO QCD for the top and bottom quark loops. 



The NLO electroweak corrections 187 



are included through the coefficient n ew and by 
neglecting terms of 0[(a 2 /A.ii){y 2 / f 2 )}. The above formula thus includes the leading 0{y 2 / f 2 ) 
corrections, as well as mixed 0[(a s /47r) 5 (a 2 /4:7r)] ones. Indeed, we assume factorisation of the 
QCD and EW corrections. Since QCD corrections are dominated by soft gluon radiation, in 



which QCD and EW effects completely factorize, this is a good approximation^ It neglects 
terms of 0[(a 2 /4vr) 2 ] and 0{v 4 /f A ). 



29 Bottom loops contribute 0(10%) to the SM decay rate and are well approximated by an effective coupling 
at the 10%-level thus leading to negligible non-factorizing contributions at the per-cent level. 
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E.3 Decay into photons 

In the SM the decay of the Higgs boson into a pair of photons is mediated by W and 



heavy fermion loops. According to the chiral Lagrangian (3.46), these two contributions 



to the total amplitude are rescaled, respectively, by the parameters Cy and c^. Similarly 
to h — > gg, the contact interaction proportional to c 77 can also contribute significantly. 
With /Air)(g 2 v 2 /M 2 ), the contribution due to the effective interaction becomes 

comparable to the loop induced contributions if {g 2 v 2 /M 2 ) 1. The partial width for a 
Higgs boson decaying into two photons implemented in eHDECAY in the framework of the 
non-linear Lagrangian is thus given by 



r(77) 



G FOi 2 em m\ 



NL 



128^ 



7P 



q=t,b,c 



+ c v A 1 (t w ) + 



47T 



-77 



(E.116) 



which is approximate at leading order in the derivative expansion, i.e. it neglects terms of 
0(m 2 h / M 2 ). Note that a em is the electromagnetic coupling in the Thomson limit, in order 
to resum large logarithms for on-shell photons. We have defined = 4m 2 /m 2 h (i = q, r, W), 
and the form factor 

A 1 (r) = - [2 + 3r + 3r (2 - r) / (r)] , (E.117) 



while Aij2{r) is given in Eq. (E.110). The top, bottom and charm quark loops receive NLO 



QCD corrections, while the effective contact interaction does not. The NLO QCD corrected 



quark form factor is denoted in Eq. (E.116) by 



A^ L 2 {r q ) = A l/2 {r q ){l + K QCD ) 



(E.118) 



where kq C d encodes the 0(a s /An) QCD corrections 89 



In the case of the SILH parametrization, the EW corrections have been incorporated as 
well. It is useful to define the SM amplitude at LO and NLO QCD level as 

^! M (77) = E t 3 ®l ^ + \<£ A ^/ 2 ^ + Al ^ ' X = LO, NLO , (E.119) 

q=t,b,c 
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where A s L q = A SM , and the deviation from the SM amplitude as 

4 /Cff \ _ „ . «m , x (C~H \ 4 



AA( 77 ) = - E I ( f + 3 « 4$° fa) - (" 



+ Cr ) 7T<2r ^1/2 fa) 



2 / 3 

The decay width implemented in eHDECAY in the SILH case is thus the following 



q=t,b,c 

2 



r(77) 



G F a 2 em m 3 



SILH 




2Re(^( 77 )^( 77 ) 



SM, 



a SM * 
NLO 



(77) ( AA( 77 ) 



327T sin 2 #jy c 7 



(E.120) 



(E.12i; 



where ^4f^ f ( 77 ) denotes the SM amplitude which comprises the NLO electroweak correc- 
tions [87l[90] . Equation (E.121) includes the leading 0(v 2 / f 2 ) and 0(m 2 /M 2 ) corrections, 



while it neglects terms of order 0(v 4 /f A ). The electroweak corrections are implemented up 
to NLO, neglecting corrections of 0[(a 2 /An)(v 2 /f 2 )} and of 0[(a 2 / An) 2 }. Finally, the QCD 
corrections are included up to NLO, and mixed terms of 0[(a 2 /An)(a s /Aii)} are neglected. 



E.4 Decay into Z 7 

In the SM the Higgs boson decay into a Z boson and a photon is mediated by a W boson 
and heavy fermion loops. Adopting the parametrization of the non-linear Lagrangian, the 
correction from the effective interaction due to the coupling c^ 7 has to be considered, too, 
and it can become as important as the loop contributions for (g 2 v 2 /M 2 ) 1. The decay 
width is therefore given by (here also a em = a em (0)): 



r(z 7 ) 



_ G 2 pOL em Tn^wi^ l 

NL ~ ^4 



m. 



x 



E 



cos 6 
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Afj 2 (r^, A^) + c v Al 1 (t w , \w) 



Air 
y/a em a 2 



c z ~ 



(E.122) 
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with Tj = 4m-/m|, Aj = 4m 2 /m|, = 21^ — sin 2 9w {i> = t,b,c,r) and the form 
factors are defined by 



Af 7 (r, A) = cos ^{4 (3 - tan 2 ^^(r, A) 



(E.123) 



+ 



1 + - ) tan 2 6 W - ( 5 4 



The functions I\ and -Z2 can be cast into the form 



h (r, A) 
h (r, A) 



r\ r 2 A 2 



2(r-A) 2(r-A) 2 
rA 



[f{r)-nX)\ 



r 2 A 



(r-A)' 



/i(r,A) 



[5 (r) - 9 (A)] 



(E.124) 



2(r-A) 



If (r) ~ f (A)] 



where /(r) is defined in Eq. (E.lll) and g(r) reads 

1 



9(r) 



\Jt — 1 arcsin 



In 



i + yr 



i-vt 



/7T 



r 



r > 1 



r < 1 



(E.125) 



The QCD radiative corrections 91 are small and thus have been neglected, while the NLO 
EW corrections are unknown. Finally, Eq. (E.122) neglects terms of 0(m 2 /M 2 ), which are 
of higher-order in the derivative expansion. 

In the SILH parametrization the decay width is computed by eHDECAY according to the 
formula 



r(z 7 ) 



SILH 



64^ 



m 



mi 



2 \ 3 



x 



I A SM (Zi) 1 2 + 2 Re ( A SM * (Zy) AA(Zj)) 



(E.126) 



+ 2 Re 



47r tan 9 



w 



chb - c H w ~ 8c 7 sin 2 6 W ) A SM *(Z^) 



y/0> em O> 2 

where we have defined the LO SM amplitude 

N C Q^ 



a sm (zj) = y: 



cosOw x l 2 



^1/2 ( r ^' A^) + ^f 7 ( T w, Aw) 



(E.127) 



55 



and the deviation from the SM amplitude 



AA( Zl ) = - ( °f + 5*) A V 2 ( T *> A*) - ( C f - 2c w ) A? (r w , X w ) . (E.128) 

\ A / COS u\y ^ ^ ' 



Equation (E.126) includes corrections of 0(v 2 / f 2 ) and 0(m 2 h / M 2 ). The EW corrections are 
unknown, and small QCD radiative corrections have been neglected. 

E.5 Decays into WW and ZZ bosons 

The Higgs boson decay into a pair of massive vector bosons is important not only above the 
threshold, but also below. For example, in the SM with irih = 125 GeV the branching ratio of 
h — > WW is about 20%. In HDECAY various options are present to compute the partial decay 
widths with on-shell or off-shell bosons, controlled by the ON-SH-WZ input parameter. In 
eHDECAY we have implemented the case ON-SH-WZ=0, which includes the double off-shell 
decays h — > W*W* , Z* Z* . For this case, which is obviously the most complete as it takes 
into account both on-shell and off-shell contributions, the partial decay width h — > V*V* 
(V = W, Z) can be written in the following compact form: 

"{m h -Qi) 2 



T(V*V*) 



dQ{ myT\ 



Jo (Ql ~ m 2 vY + m 2 v T 2 v 



dQ 2 , myY\ 



(Q 2 2 - m 2 v Y + m 2 v T 2 v 



T(VV), (E.129) 



where Q 2 , Q\ are the squared invariant masses of the virtual gauge bosons and my and T\ 



Cbyy — C VV 



mt 
m 



2 ' 

V 



their masses and total decay widths. In the parametrization of Eq. (3.46), by defining 

_ c v9v m\ f-Ri^m 
avav - ~y~ ' (E.ldO) 

the squared matrix element r(Kl^) reads 
r(W)L r = T SM (VV) x { c 2 v - 2r v 



Cbyy 



1 - 



Q\ + Ql\ Ql + Ql 



mt 



+ &VdV~ 



mt 



+ Cyayy 



XjQlQlm 2 ) (l-(Cft + QD/m l) 
A (Ql, Q 2 , m 2 h ) + 12 Q\Q\j m\ 



(E.i3i; 
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T SM (VV) = 5 -^^\lHQlQlrn 2 h ) (A (Ql Ql m\) + 12 ^j , (E.132) 
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where 5y = 2(1) for V = W(Z) and X(x,y,z) = (1 — x/z — y/z) 2 — Axy/z 2 . The second 



and third term in Eq. (E.131) represent the interference between the tree-level contribution 



and the one from the derivative operators. They are of order 0(m\/M 2 ), hence next-to- 
leading in the chiral expansion compared to the tree-level contribution; we have consistently 
neglected terms quadratic in ayv and avav, since they are of 0(m^/M 4 ), which is beyond 



the accuracy of the effective Lagrangian (3.46) 



In the SILH parametrization the squared matrix element T(l A V A ) implemented in eHDECAY 
reads 



Y{VV) 



^SILH 



(VV) + T SM {VV) 



SILH 



Re (Af M A s J) 



(E.133) 



where Aq M denotes the SM LO amplitude and A^ 1 is the SM amplitude which comprises 



the NLO EW corrections [93] (the same remark as in footnote 21 applies). Furthermore, 

Qi + Ql\ , - Ql + QY 



^SILH, 



(VV) = T SM (VV) x I 1 - c H - 2 
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mi 
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(E.134) 



(E.135) 
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a Z z = —2 — y (chw + chb tan 2 9w — 4c 7 tan 2 9w sin 2 ^) 



&zdz — —2 



m h 
2m| 



(c w + c HW + (c B + c HB ) tan 2 6>iy) 



(E.136) 



The amplitude (E.133) includes terms of 0(v 2 /f 2 ), 0(m 2 h / M 2 ) and 0(«2/47r), while it ne- 
glects contributions of 0(t> 4 // 4 ) and 0[(a2/47r) 2 ]. As explained in Section 



4.2 



the correc- 



tions of 0[(a;2/47r)(t> 2 // 2 )] are only partly included through the terms proportional to Chw 
and chb- 
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E.6 Numerical formulas for the decay rates in the SILH Lagrangian 

We display here numerically approximated formulas of the Higgs decay rates valid at linear 



order in the effective coefficients cij of the SILH Lagrangian (2.2) for rrih = 125 GeV. Although 
these numerical results have been partly shown in Section |3j we think it is useful to collect 
them here in a systematic and complete way All the ratios T/Tsm have been computed by 
switching off the EW corrections, since their effect on the numerical pref actor appearing in 
front of the coefficients q is of order (v 2 / f 2 )(ct2/4:7r) and thus beyond the accuracy of the 
formulas implemented in eHDECAY. Conversely, we have fully included the QCD corrections, 
as they multiply both the SM and the NP terms. The numerical results are thus the following: 



l-c H -2c0, (E.137) 

l-c H + 2.2c w + 3.7c HW , (E.138) 

, - 1 - c H + 2.0 (c w + tan 2 9 w c B ) 

T{h ->■ Z(*)Z*)sm (E.139) 

+ 3.0 (chw + tan 2 #vy chb) — 0.26 c 7 , 



r(h ->■ w^w*) 

T(h -+ W^W*) SM 



T(h -» Zi) 



r(h 77) 



1 - c H + 0.12 q - 5 ■ 10" 4 c c - 0.003 c b - 9 ■ 10~ 5 c T 
+ 4.2 cV + 0.19 (chw — chb + 8 c 7 sin 2 ^) 



(E.140) 



em 



r(/i -> 77 ) 
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1 - c H + 0.54 q - 0.003 c c - 0.007 c b - 0.007 c T 

47T 

+ 5.04c w - 0.54 c. 



(E.141) 
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-»■ qq) Att 

' -1 -ch- 2.12 c 4 + 0.024 c c + 0.1 c b + 22.2 c g — . (E.142) 



T(h^gg) SM ' ' ' ' " « 2 
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F How to run eHDECAY: Input /Output Files 

The program eHDECAY is self-contained, like the original code HDECAY on which it is based. 
All the new features related to the Lagrangian parametrizations proposed in this paper are 
encoded in the main source file, ehdecay.f, while other linked routines are taken over from the 
original version. Of course eHDECAY, besides calculating Higgs branching ratios and decay 
widths according to the non-linear, SILH or MCHM4/5 Lagrangian, also calculates the SM 
and MSSM ones, exactly as HDECAY 5.10 does. The choice can be done through the flags 
HIGGS and COUPVAR set in the input file. The input file for eHDECAY has been called 
ehdecay.in and is based on the file hdecay.in of the official version 5.10, supplemented by 
further input values. Thus, with the flag LAGPARAM the user can choose between the 



general SILH parametrization Eq. (2.2), the model-specific parametrizations MCHM4 and 



MCHM5 and the general non- linear Lagrangian parametrization Eq. (3.46). Furthermore 



the various related couplings can be set. The input values are explained in the following: 

COUPVAR, HIGGS : If HIGGS=0 and COUPVAR=l, then the Higgs decay widths and 
branching ratios are calculated within the parametrization chosen by: 

LAGPARAM: 



Non-linear Lagrangian parametrization Eq. (3.46) 



SILH parametrization Eq. (2.2) 



MCHM4/5 parametrization (cf. Table Q 



IELW : Turn off (0) or on (1) the electroweak corrections for the SILH, MCHM4 or MCHM5 
parametrization. Fj 

For the non-linear Lagrangian the following parameters have to be set for the couplings of 
the various vertices: 



31 



30 Note, that this parameter IELW has nothing to do with the parameter ELWK in the input file of 
HDECAY, where the meaning of this flag is different. 
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We explain them here all, although they are in part already present in the input file for HDECAY 5.10. 
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CV : hVV vertex, (V=W, Z) Ctau: hrr vertex 



Cmu : hjiji vertex 



Ct: hit vertex 



Cb: hbb vertex 



Cc: hcc vertex 



Cs: hss vertex 



Cgaga: coupling c, 



•77 



Cgg: coupling c, 



CZga: coupling c z . 



7 



CWW : coupling Cy/w CZZ : coupling c zz 



CWdW : coupling c W g W 



CZdZ : coupling c^g^ 



In case of the SILH parametrization the input values to be set in order to calculate the 
various couplings are: 

CHbar : ch Ctaubar : c T Cmubar : c M Ctbar : c t Cbbar : q, 



CHBbar : chb Cgambar: c 7 Cgbar: c g 
In the MCHM4/5 parametrization we have the input values: 

FERMREPR : 
1: MCHM4 
2: MCHM5 

XI : the value for £ 

For example: 

COUPVAR = 1 
HIGGS = 

************** LAGRANGIAN - chiral 1 - SILH 2 - MCHM4/5 ************** 
LAGPARAM = 

**** Turn off (0) or on (1) the elw corrections for LAGPARAM = 1 or 2 **** 



Ccbar: c c 



Csbar : c s 



CWbar : cw CBbar : cb CHWbar : chw 
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IELW = 1 

******************* VARIATION OF HIGGS COUPLINGS************************* 

CW = l.DO 

CZ = l.DO 

Ctau = 0.95D0 

Cmu = 0.95D0 

Ct = 0.95D0 

Cb = 0.95D0 

Cc = 0.95D0 

Cs = 0.95D0 

Cgaga = 0.005D0 

Cgg = O.OOIDO 

CZga = 0.D0 

CWW = O.DO 

CZZ = O.DO 

CWdW = O.DO 

CZdZ = O.DO 

computes the branching ratios for cy = 1, = 0.95 (ip = t, b, c, s, r, fi), c 77 = 0.005, 
Cgg = 0.001 and c^ 7 = cww = czz = cwaw = czaz = in the general parametrization 



Eq. (3.46). The output is written into the files br.effl and br.eff2, where the Higgs mass, 
branching ratios and total width are reported. For the previous example, at = 125 GeV 
and for all the other parameters set at their standard values, the output reads 

MHSM BB TAU TAU MU MU SS CC TT 



125.000 0.5895 0.5654E-01 0.2002E-03 0.2161E-03 0.2569E-01 0.000 



MHSM GG GAM GAM Z GAM WW ZZ WIDTH 
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125.000 0.9611E-01 0.1932E-03 0.1526E-02 0.2045 0.2554E-01 0.4129E-02 



All the input parameters of the corresponding run are printed out in the file br. input. Other- 
wise, setting COUPVAR=0, the program produces the usual output files with SM or MSSM 
results according to the HDECAY 5.10 version. 
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